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Abstract 

We study the decay constants and form factors of the ground-state s-wave and low-lying p-wave 
mesons within a covariant light-front approach. Numerical results of the form factors for transi- 
tions between a heavy pseudoscalar meson and an s-wave or p-wave meson and their momentum 
dependence are presented in detail. In particular, form factors for heavy-to-light and B — > D** 
transitions, where D** denotes generically a p- wave charmed meson, are compared with other model 
calculations. The experimental measurements of the decays B~ D**tt~ and B DD** are em- 
ployed to test the decay constants of D** and the B D** transition form factors. The heavy 
quark limit behavior of the decay constants and form factors is examined and it is found that the 
requirement of heavy quark symmetry is satisfied. The universal Isgur-Wise (IW) functions, one 
for s-wave to s-wave and two for s-wave to p-wave transitions, are obtained. The values of IW 
functions at zero recoil and their slope parameters can be used to test the Bjorken and Uraltsev 
sum rules. 
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I. INTRODUCTION 



Mesonic weak transition form factors and decay constants are two of the most important ingredi- 
ents in the study of hadronic weak decays of mesons. There exist many different model calculations. 
The light-front quark model P, |^ is the only relativistic quark model in which a consistent and fully 
relativistic treatment of quark spins and the center-of-mass motion can be carried out. This model 
has many advantages. For example, the light-front wave function is manifestly Lorentz invariant as 
it is expressed in terms of the momentum fraction variables in analog to the parton distributions in 
the infinite momentum frame. Moreover, hadron spin can also be correctly constructed using the 
so-called Melosh rotation. This model is very suitable to study hadronic form factors. Especially, 
as the recoil momentum increases (corresponding to a decreasing q'^), we have to start considering 
relativistic effects seriously. In particular, at the maximum recoil point = where the final-state 
meson could be highly relativistic, there is no reason to expect that the non-relativistic quark model 
is still applicable. 

The relativistic quark model in the light-front approach has been employed to obtain decay 
constants and weak form factors 0,0,0,0,0] • There exist, however, some ambiguities and even some 
inconsistencies in extracting the physical quantities. In the light-front quark model formulation 
one often picks up a specific Lorentz frame (e.g. the purely longitudinal frame q± = 0, or the 
purely transverse frame q^ = q^ + q"^ = 0) and then calculates a particular component (the "plus" 
component) of the associated current matrix element. Due to the lack of relativistic covariance, 
the results may not be unique and may even cause some inconsistencies. For example, it has 
been pointed out in 0| that in the q± = frame, the so-called Z-diagram contributions must be 
incorporated in the form-factor calculations in order to maintain covariance. Another issue is that 
the usual recipe of taking only the plus component of the current matrix elements will miss the 
zero- mode contributions and render the matrix element non-covariant. A well known example is the 
electromagnetic form factor i*2('7^) of the vector meson (see e.g. j^). In other words, the familiar 
expression of /y, for example, in the conventional light-front approach ^ is not trustworthy due to 
the lack of the zero-mode contributions. As a consequence, it is desirable to construct a covariant 
light-front model that can provide a systematical way of exploring the zero-mode effects. Such a 
covariant model has been constructed in 0] for heavy mesons within the framework of heavy quark 
effective theory. 

Without appealing to the heavy quark limit, a covariant approach of the light-front model for the 
usual pseudoscalar and vector mesons has been put forward by Jaus (for a different approach, 
see lul). The starting point of the covariant approach is to consider the corresponding covariant 
Feynman amplitudes in meson transitions. Then one can pass to the light-front approach by using 
the light-front decomposition of the internal momentum in covariant Feynman momentum loop 



integrals and integrating out the p~ = — component [1^. At this stage one can then apply 
some well-studied vertex functions in the conventional light-front approach after p~ integration. 
It is pointed out by Jaus that in going from the manifestly covariant Feynman integral to the 
light-front one, the latter is no longer covariant as it receives additional spurious contributions 
proportional to the lightlike vector u;^ = (1,0,0,-1). This spurious contribution is cancelled after 
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correctly performing the integration, namely, by the inclusion of the zero mode contribution so 
that the result is guaranteed to be covariant. Before proceeding, it is worth mentioning that in the 
literature there is a controversy about the zero mode contributions to the vector decay constant fy 
and the form factor Ai{q^) in the pseudoscalar to vector transition: While Jaus llOl. llJ claimed that 
there are zero effects in the aforementioned two quantities, Bakker, Choi and Ji [ajlUl argued that 
both fv and Ai^q"^) are free of zero-mode contributions. This issue will be addressed in Sec. III.B. 

The main purposes of this work are twofold: First, we wish to extend the covariant analysis of 
the light-front model in 1^ to even-parity, p-wave mesons. Second, the momentum dependence of 
the form factors is parametrized in a simple three-parameter form so that the reader is ready to use 
our numerical results as the analytic expressions of various form factors in the covariant light-front 
model are usually complicated (see Sec. III). Interest in even-parity charmed mesons has been 
revived by recent discoveries of two narrow resonances: the 0+ state Z?*q(2317) 0] and the Pi^'^ 



state L>si(2460) and two broad resonances, Z?j5(2308) and Di(2427) IjJ.^ Furthermore, the 
hadronic B decays such as B ^ D**tt and B D**D have been recently observed, where D** 
denotes a p-wave charmed meson. A theoretical study of them requires the information of the 
B D** form factors and the decay constants of D** and D** . In the meantime, three body 
decays of B mesons have been recently studied at the B factories: BaBar and Belle. The Dalitz 
plot analysis allows one to see the structure of exclusive quasi-two-body intermediate states in 
the three-body signals. The p-wave resonances observed in three-body decays begin to emerge. 
Theoretically, the Isgur-Scora-Grinstein-Wise (ISGW) quark model [l9| is so far the only model 
in the literature that can provide a systematical estimate of the transition of a ground-state s- 
wave meson to a low-lying p-wave meson. However, this model and, in fact, many other models in 
P ^ P,V {P: pseudoscalar meson, V: vector meson) calculations, are based on the non-relativistic 
constituent quark picture. As noted in passing, the final-state meson at the maximum recoil point 
= ov in heavy-to-light transitions could be highly relativistic. It is thus important to consider 
a relativistic approach. 

It has been realized that the zero mode contributions can be interpreted as residues of virtual pair 
creation processes in the g+ — > limit the calculation of the zero mode contribution 

is obtained in a frame where the momentum transfer q~^ vanishes. Because of this (g"*" = 0) 
condition, form factors are known only for spacelike momentum transfer q"^ = —qj_ < 0. One needs 
to analytically continue them to the timelike region 0|, where the physical decay processes are 
relevant. Recently, it has been shown that within a specific model, form factors obtained directly 
from the timelike region (with q^ > 0) are identical to those obtained by the analytic continuation 
from the spacelike region 

There are some theoretical constraints implied by heavy quark symmetry (HQS) in the case 
of heavy-to- heavy transitions and heavy-to- vacuum decays [2]|. It is important to check if the 
calculated form factors and decay constants do satisfy these constraints. Furthermore, under HQS 
the number of the independent form factors is reduced and they are related to some universal 



^ We follow the naming scheme of the Particle Data Group [l^l to add a superscript "*" to the states if the 
spin-parity is in the "normal" sense, J-^ = 0+, 1^, 2+, • • •. 
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FIG. 1: Feynman diagrams for (a) meson decay and (b) meson transition amplitudes, where P'(") 
is the incoming (outgoing) meson momentum, p^"^ is the quark momentum, p2 is the anti-quark 
momentum and X denotes the corresponding V — A current vertex. 



Isgur-Wise (IW) functions. In this work, we shah follow [9| to evaluate the form factors and decay 
constants in a covariant light-front formulism within the framework of heavy quark effective theory. 
It is found that the resultant decay constants and form factors do agree with those obtained from 
the covariant light-front approach and then extended to the heavy quark limit. The relevant IW 
functions, namely, ^, t^i2 are obtained. One can then study some properties of these IW 



functions, including the slopes and sum rules 22, 



The paper is organized as follows. In Sec. II, we give the calculations for the decay constants 
of s-wave and p-wave mesons in a covariant light-front model. The calculation for s-wave meson 
transitions has been done by Jaus 01 • We extend it to the p-wave meson case. In Sec. Ill, P — > 
P ,V ,A,T {S,A,T standing for scalar, axial-vector and tensor mesons, respectively) transitions are 
considered. It is interesting to notice that the analytic forms of P — > S, A transitions are similar to 
that of P ^ P,V transitions, respectively, while the P ^ T calculation needs formulas beyond [lo| . 
We provide numerical results for B and D decay form factors and their q"^ dependence. These results 
are then compared to the other model calculations. In Sec. IV, properties of the decay constants 
and form factors in the heavy quark limit are studied. The universal Isgur-Wise functions, one 
for s-wave to s-wave and two for s-wave to p-wave transitions, are obtained. Their values at zero 
recoa and their dope parameters can be used to test the sun, rules derrved by Bjorken Q and by 
Uraltsev [2^. Conclusion is given in Sec V followed by two Appendices devoted to the derivations 
of conventional light-front vertex functions and some useful formulas. 

II. FORMALISM OF A COVARIANT LIGHT-FRONT MODEL 
A. Formalism 

In the conventional light-front framework, the constituent quarks of the meson are required to 
be on their mass shells (see Appendix A for an introduction) and various physical quantities are 
extracted from the plus component of the corresponding current matrix elements. However, this 
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procedure will miss the zero- mode effects and render the matrix elements non-covariant. Jaus [lO|] 
has proposed a covariant light-front approach that permits a systematical way of dealing with the 
zero mode contributions. Physical quantities such as the decay constants and form factors can be 
calculated in terms of Feynman momentum loop integrals which are manifestly covariant. This 
of course means that the constituent quarks of the bound state are off-shell. In principle, this 
covariant approach will be useful if the vertex functions can be determined by solving the QCD 
bound state equation. In practice, we would have to be contended with the phenomenological 
vertex functions such as those employed in the conventional light-front model. Therefore, using the 
light-front decomposition of the Feynman loop momentum, say p^, and integrating out the minus 
component of the loop momentum p~, one goes from the covariant calculation to the light-front 
one. Moreover, the antiquark is forced to be on its mass shell after p~ integration. Consequently, 
one can replace the covariant vertex functions by the phenomenological light-front ones. 

As stated in passing, in going from the manifestly covariant Feynman integral to the light-front 
one, the latter is no longer covariant as it can receive additional spurious contributions proportional 
to the lightlike four vector ui. The undesired spurious contributions can be eliminated by the 
inclusion of the zero mode contribution which amounts to performing the p~~ integration in a 
proper way in this approach. The advantage of this covariant light-front framework is that it 
allows a systematical way of handling the zero mode contributions and hence permits to obtain 
covariant matrix elements. 

To begin with, we consider decay and transition amplitudes given by one-loop diagrams as shown 
in Fig. ^ for the decay constants and form factors of ground-state s-wave mesons and low- lying p- 
wave mesons. We follow the approach of and use the same notation. The incoming (outgoing) 
meson has the momentum = p'^^"^ -\- p2 , where p'j"^ and p2 are the momenta of the off-shell 
quark and antiquark, respectively, with masses m'/"^ and m2- These momenta can be expressed in 
terms of the internal variables {xi,p'j_), 

Pi^2 = ^i,2P'^ , Pl,2± = xi,2P±±P±, (2.1) 
with XI + X2 = 1. Note that we use P' = (P'+ , P'" , P[ ) , where P'^ = ± P'^, so that P'^ = 
P'^P'^ — P'^. In the covariant light-front approach, total four momentum is conserved at each 
vertex where quarks and antiquarks are off-shell. These differ from the conventional light-front 
approach (see, for example 0,01) where the plus and transverse components of momentum are 
conserved, and quarks as well as antiquarks are on-shell. It is useful to define some internal 
quantities analogous to those defined in Appendix A for on-shell quarks: 



Xl X2 V " 



- ym^ +P^+P,, Pz- 2 2x2M^ ■ ^ ' 

Here can be interpreted as the kinetic invariant mass squared of the incoming qq system, and 

ei the energy of the quark i. 

It has been shown in [l^ that one can pass to the light-front approach by integrating out the 
component of the internal momentum in covariant Feynman momentum loop integrals. We 



5 



need Feynman rules for the meson-quark-antiquark vertices to calculate the amplitudes shown in 
Fig. 1. These Feynman rules for vertices (i^'m) of ground-state s-wave mesons and low- lying p-wave 
mesons are summarized in Tabled As we shall see later, the integration of the minus component 
of the internal momentum in Fig. 1 will force the antiquark to be on its mass shell. The specific 
form of the covariant vertex functions for on-shell quarks can be determined by comparing to the 
conventional vertex functions as shown in Appendix A. Next, we shall use the decay constants as 
an example to illustrate a typical calculation in the covariant light-front approach. 

B. Decay constants 

The decay constants for J = 0, 1 mesons are defined by the matrix elements 

{0\A^\P{P')) ^ < = ifpP',, {0\V^\S{P')) ^Af, = fsP',, (2.3) 

where the ^'^"'"^Lj = ^5o, '^Pq, ^"S*!, ^Pi, ^Pi and ^P2 states of (^iq2 mesons are denoted by P, 5, V , 
^A, and T, respectively. Note that a ^P2 state cannot be produced by a current. It is useful to 
note that in the SU(N)-flavor limit {m'^ = 1712) we should have vanishing fs and /i^. The former 
can be seen by applying equations of motion to the matrix element of the scalar resonance in Eq. 
((T^ to obtain 

mlfs = i{m[ - m2){0\qiq2\S). (2.4) 

The latter is based on the argument that the light ^Pi and ^Pi states transfer under charge conju- 
gation as 

M,^(3pi) ^ M,"(3pi), M'XPi)^-M^CPi), (a = 1,2,3), (2.5) 

where the light axial-vector mesons are represented by a 3 x 3 matrix. Since the weak axial- 
vector current transfers as (A/^)^ {^fj.)b under charge conjugation, it is clear that the decay 

TABLE I: Feynman rules for the vertices (iF^^) of the incoming mesons-quark-antiquark, where 
p'l and p2 are the quark and antiquark momenta, respectively. Under the contour integrals to be 
discussed below, H'j^j and W^j are reduced to h'f^j and w'j^^, respectively, whose expressions are 
given by Eq. (|2.11|) . Note that for outgoing mesons, we shall use i(7or'2f7o) for the corresponding 
vertices. 







pseudoscalar (^Sq) 
vector {"^Si) 

scalar (^Pq) 
axial (3Pi) 

axial (iPi) 

tensor (^P2) 


iH'yl-if, - i^{p'i -P2)^^ 
-iH's 

-iHi^hf, + wriPi -P2)^^h5 
i^H^bii - w^^ip'i - P2)^,]{p'l - P2)u 
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constant of the ^Pi meson vanishes in the SU(3) hmit 24 1. This argument can be generahzed to 
heavy axial-vector mesons. In fact, under similar charge conjugation argument [(V"^)a — > —{Vfj.)b^ 
M^(^Po) ~^ ^bi^^o)] one can also prove the vanishing of fs in the SU(N) limit. 

Furthermore, in the heavy quark limit (m'^ — > oo), the heavy quark spin sq decouples from the 
other degrees of freedom so that sq and the total angular momentum of the light antiquark j are 
separately good quantum numbers. Hence, it is more convenient to use the L-^j = P2^'^ , Pi^"^ : Pi^'^ 



and Pn^^ basis. It is obvious that the first and the last of these states are and ^Pq, respectively, 
while H 



P 



,3/2 



Pl) + 



1 

7! 



'Pi 



p 



1/2 



1 

7! 



Heavy quark symmetry (HQS) requires (see Sec. IV) 
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'Pi 



'Pi 



(2.6) 



fv = fp, /ai/2 = fs, /a3/2 = 0, (2.7) 

where we have denoted the Pi^"^ and Pj^^^ states by A^^'^ and A^^'^, respectively. These relations in 
the above equation can be understood from the fact that (5^^, 5^^^), (Pq^'^^, Pj'^'^^) and (Pf''^, P2 ''^) 
form three doublets in the HQ limit and that the tensor meson cannot be induced from the V — A 
current. It is important to check if the calculated decay constants satisfy the non-trivial SU(N)- 
flavor and HQS relations. 

We now follow [lfi| to evaluate meson decay constants. The matrix element for the annihilation 
of a pseudoscalar state via axial currents can be easily written down and it has the expression 

.P 2 r A, H'p p 

{2tt)^J ^'N{N2 

where 



(2i 



•s^ = Tr[7M75(/i + m'i)75(- ^2 + m2)] 

= -A[m[P'^ + {m2-m[)p[^], (2.9) 

N[ = p'l — m'l + ie and A^2 = P2 ~ "'-i + We need to integrate out p'{' in A^. As stressed in 0], 
if it is assumed that the vertex function H' has no pole in the upper complex p'l' plane, then the 
covariant calculation of meson properties and the calculation of the light-front formulism will give 
identical results at the one-loop level. Therefore, by closing the contour in the upper complex 
plane and assuming that H'p is analytic within the contour, the integration picks up a residue at 
P2 = P2, where = The other momentum is given by momentum conservation, p'l = P' — p2. 
Consequently, one has the following replacements: 

n'l = xi{M'^- 

M — H'm{Pi,P2) = ^M; 
W'Mip'lp^)^^'^^, 



N'l N'l 



-12 
Pi 



M, 



i2\ 



H', 



M 



N[N2 



H 



M 



-ITT 



dx2d Pj_ 

X2N[ 



h', 



(2.10) 
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in a generic one- loop vacuum to particle M amplitude . In this work the explicit forms of h'j^ 
and w'^j are given by (see Appendix A) 

12 nrl2\ /^i^2 1 _/ 



h'p = h'y = (M'2 - M'i)J-^-^^' 

h's = ./|/i'3^ = (M'2-M^2) / 



wy = M'q + m\ + ma, w'zj^ = — ^ , w[, = 2 , (2.11) 

m,]^ — m2 

where 99' and 99^ are the light-front momentum distribution amplitudes for s-wave and p-wave 
mesons, respectively. There are several popular phenomenological light-front wave functions that 
have been employed to describe various hadronic structures in the literature. In the present work, 
we shall use the Gaussian-type wave function '27] 



93 = {X2,p^) 




p'z+p'T 



2/3' 



12 



c,; = ^;(X2,A) = ^ = (2.12) 

The parameter (3' is expected to be of order Aqcd- The derivation of these vertex functions is 
shown in Appendix A. 

The matrix element A^ can be evaluated readily by using above equations. However, A^ ob- 
tained in this way contains a spurious contribution proportional to uj^ = (a)~,u;+,a)_L) = (2,0, 0_l). 
It arises from the momentum decomposition of pf^ 

n'r = ip' - 



Pl = -P2 

X 



In fact, after the integration, p'^ can be expressed in terms of two external vectors, P' and a). 

\M 

' - ^''\p:, + ^.^Ap'-p\-^p" 



Therefore, in the integrand of one has 



= ^iK + iT^^'^f^[-N2 + N[+ mf -ml + {l- 2xi)M'^]. (2.14) 

The symbol = in the above equation reminds us that it is true only in the A^ integration. There 
is one missing piece in the contour integration, namely, the contribution of the zero mode from the 
p^ = region The appearance of N2 in the numerator as shown in the above equation (|2.8|) 
also prompts an extra care in performing the p'{' contour integration. It is interesting that this 
zero mode contribution provides a cue for the spurious term in Aff . As shown in ^^1) the inclusion 
of the zero mode contribution in A^^ matrix elements in practice amounts to the replacements 

p[ xiP', N2^N[+ mf -ml + {l- 2xi)M'^, (2.15) 
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in the under the integration. By virtue of Eqs. (|2.3j) . (|2.1U|) and 1)2. 15j) . we obtain 0| 



167r3 



h' 



4(m'^X2 + m2Xi). 



(2.16) 



It should be stressed that fp itself is free of zero mode contributions as its derivation does not involve 



the replacement of N2 (see also Sec. III.B). With the explicit form of h'p shown in Eq. ()2.1H) . the 



familiar expression of fp in the conventional light-front approach 



, namely, 



ff 



' 167r3 



dx2(f'p_i_- 



1 



jj {m[x2 + m2Xi) ip'{x2,p±), 



(2.17) 



^XiX2M(^ 

is reproduced. 

The decay constant of a scalar meson can be obtained in a similar manner. By using the 
corresponding Feynman rules shown in Tabled we have 



At 



I ^Vi^T¥[7^(/i + m[){-i){- h + "^2)]. (2.18) 



Note that the trace (= s^) in the above equation is related to in Eq. 1)2. 8|) by the replacement of 
m2 — > —m2 and by adding an overall factor of —i. Likewise, by using Eqs. (|2.3() . 1)2.10(1 and (|2.15() . 
it follows that 



h'o 



-A{m'iX2 — m2Xi). 



(2.19) 



iX2[M'^-M(^) 

For m']^ = 777-2, the meson wave function is symmetric with respect to xi and 3:2, and hence = 0, 
as it should be. 

We now turn to the decay constants of vector and axial- vector mesons. The decay amplitude 
for a vector meson is given by 

{p'l - P2) 



A 



V 



.2 

' (2vr)^ 



We consider the case with the transverse polarization 

/ 2 



e(±) 



1 



(- ^2 + m2) ) e ' 



(2.20) 



(2.21) 



Contracting AY, with e*{±) and applying Eqs. (ITTn|) and TTT^ lead to 0^ 



Iv 



dx2d'^p'± 



h' 



V 



a;iX2(M'2 - M^2^ 



Mil III \ 12 , 

^ Q — mi[mi — 1712) — p _\_ + 



,2 , m'l + m2 



12 
■P± 



(2.22) 



^ When is contracted with the longitudinal polarization vector e^(0), fv will receive additional contri- 
butions characterized by the B functions defined in Appendix B (see Eq. (3.5) of 14J) which give about 
10% corrections to fv for the vertex function h'y used in Eq. (|2.11() . It is not clear to us why the result 
of fv depends on the polarization vector. Note that the new residual contributions are absent in the 
approach of in which a different scheme has been developed to identify the zero mode contributions 
to the decay constants and form factors. 
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We wish to stress that the vector decay constant obtained in the conventional hght-front model 
does not coincide with the above result (|2.22|1 owing to the missing zero mode contribution, whose 
presence is evidenced by its involvement of N2 !3) Qj- Since Jl^ {^i^) is related to by a 
suitable replacement of H'y —H'^^^i^^ and m2 —Tn2, Wy — > — M^3^(i^) in the trace (only the 
1/1^' terms being kept in the "^A case), this allows us to readily obtain 



j;iX2(M'2 - M^' 



/2\ 



xiMq - mi (mi + 7712) - Pi V\ 

f ^2 / ( ^'1 - ^2 ,2' 



f--I^J "-^V. j^^^P^ ^ (2.23) 



It is clear that fiy^ = for ni[ = m2. The SU(N)-flavor constraints on fs and fi^ are thus satisfied. 
The HQS relations on decay constants will be discussed in Section IV. 

In order to have a numerical study for decay constants, we need to specify the constituent quark 
masses and the parameter B appe aring in the Gaussian- type wave function H2.12() . For constituent 
quark masses we use 0, 0, liol. I29I 



m„,rf = 0.26GeV, m^ = 0.37GeV, mc = 1.40GeV, mb = 4.64GeV. (2.24) 

As we shall see in Sec. Ill, the masses of strange and charmed quarks are constrained from the 
measured form-factor ratios in semileptonic D K*ii' decays. Shown in Tables O and IIIII are 
the input parameter f3 and decay constants, respectively. In Table IIIII the decay constants in 
parentheses are used to determine /?. For the purpose of an estimation, for j?-wave mesons in D, 
Dg and B systems we shall use the p»amet„s obtained in the IS0W2 model Q. the improved 
version of the ISGW model, up to some simple scaling. Several remarks are in order: (i) The 
values of the parameter Pv presented in Table ^ are slightly smaller than the ones obtained in the 
earlier literature. For example, Pp = 0.26, Pk* = 0.27 and Pd* = 0.38 are obtained here using 
the Gaussian-type wave function, while the corresponding values are 0.30,0.31,0.46 in 0]. This 
is because we have utilized the correct light-front expression for the vector decay constant fv [cf. 
Eq. ()2.22|) ]. It is interesting to notice that Pv in the ISGW2 model also has a similar reduction 
due to hyperfine interactions, which have been neglected in the original ISGW model in the mass 
spectrum calculation, (ii) The P parameters for p-wave states of D, Ds and B systems are the 



TABLE II: The input parameter P (in units of GeV) in the Gaussian- type wave function (|2.12|) . 



25+1^^ 


f^ud 


Psu 


Pcu 


Pes 


Pbu 


'So 


0.3102 


0.3864 


0.4496 


0.4945 


0.5329 


'Sr 


0.2632 


0.2727 


0.3814 


0.3932 


0.4764 


'Po 


Pa, 


Pk{^p,) 


0.3305 


0.3376 


0.4253 


'Pi 


0.2983 


0.303 


0.3305 


0.3376 


0.4253 


'Pi 


Pa, 


Pk{^Px) 


0.3305 


0.3376 


0.4253 
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TABLE III: Mesonic decay constants (in units of MeV) obtained by using Eqs. H2.16() . (|2.19() . 



(|2.22|) and (|2.23j) . Those in parentheses are taken as inputs to determine the corresponding /3's 
shown in Table ITTl The decay constant fKi{i270) — l'''^ MeV is also used as an input (see the text 
for detail). 





f ud 


fsu 


feu 


fes 


fbu 


'So 


(131) 


(160) 


(200) 


(230) 


(180) 


'S^ 


(216) 


(210) 


(220) 


(230) 


(180) 


'Po 





21 


86 


71 


112 


'Pi 


(-203) 


-186 


-127 


-121 


-123 


'Pi 





11 


45 


38 


68 


pl/2 






130 


122 


140 


p3/2 






-36 


-38 


-15 



smallest when compared to f3py. (iii) The decay constants of ^Pi and states have opposite 
signs to that of ^Pi or p'^"^ as can be easily seen from Eq. p.6|) . 

In principle, the parameter (3 for p-wave mesons can be determined from the study of the meson 
spectroscopy. Although we have not explored this issue in this work, it is important to keep in mind 
that /3's are closely related to meson masses. In Table HTll we have employed \fai \ = 203 MeV and 
foi = fos ^ inputs. It is generally argued that ai(1260) should have a similar decay constant as 
the p meson. Presumably, fai can be extracted from the decay r — > ai(1260)z^T-- Though this decay 
is not shown in the Particle Data Group (PDG) an experimental value of \fax \ = 203 it 18 MeV 
is nevertheless quoted in ,31].^ Contrary to the non-strange charmed meson case where D* has a 
slightly larger decay constant than D, the recent measurements of P ^ oi*^ D^*") 
that the decay constants of D* and Dg are similar. Hence we shall take fo* = fn„. As for the 
decay constant of P*, a recent lattice calculation yields fs* / fs = 1-01 ± O.Ollo;!]^ Therefore 
we will set fs* = fs in Table ITTll 

It is clear from Eq. ()2.4|) that the decay constant of light scalar resonances is largely suppressed 
relative to that of the pseudoscalar mesons owing to the small mass difference between the con- 
stituent quark masses. However, as shown in Table IIIII this suppression becomes less restrictive 
for heavy scalar mesons because of heavy and light quark mass imbalance. Note that what is 
the underlying quark structure of light scalar resonances is still controversial. While it has been 
widely advocated that the light scalar nonet formed by (t(600), k(800), /o(980) and ao(980) can 
be identified primarily as four-quark states, it is generally believed that the nonet states /o(1370), 
an (1450), -K'o(1430) and /o(1500)//o(1710) are the conventional qif states (for a review, see e.g. 
[3J|). Therefore, the prediction of fs = 21 MeV for the scalar meson in the su content (see Table 
IIII|) is most likely designated for the Kq{1430) state. Notice that this prediction is slightly smaller 



The decay constant of ai can be tested in the decay 5+ — + D^a^ which receives the main contribution 
from the color-allowed amplitude proportional to fa^F^^ {ra^^_^). 



18, 31 indicate 
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than the result of 42 MeV obtained in |35|] based on the finite-energy sum rules, and far less than 
the estimate of (70 it 10) MeV in It is worth remarking that even if the light scalar mesons are 
made from 4 quarks, the decay constants of the neutral scalars cj(600), /o(980) and ag(980) must 
vanish owing to charge conjugation invar iance. 

In principle, the decay constant of the scalar strange charmed meson D*q can be determined 
from the hadronic decay B — > DD*q since it proceeds only via external VF-emission. Indeed, a 
recent measurement of the DD*q production in B decays by Belle 0] indicates a fD*^ of order 60 
MeV 138'] which is close to the expectation of 71 MeV (see Table HTH) . In Sec. III.E we will discuss 
more about DD** productions in B decays. The smallness of the decay constant fD*^ relative to 
can be seen from Eqs. H2.16() and H2.19() that 

/d^d^o) « J dx2---[incX2±ms{l- X2)]. (2.25) 

Since the momentum fraction X2 of the strange quark in the Ds{D*q) meson is small, its effect 
being constructive in Dg case and destructive in D*q is sizable and explains why fD*^^/ fos ~ 0-3- 

Except for ai and bi mesons which cannot have mixing because of the opposite C-parities, 
physical strange axial-vector mesons are the mixture of ^Pi and ^Pi states, while the heavy axial- 

1/2 3/2 

vector resonances are the mixture of P^ and P^ . For example, 7^^1(1270) and i^i(1400) are the 
mixture of Kzp^ and Kip_^ (denoted by Kia and Kib, respectively, by PDG 0]) owing to the mass 
difference of the strange and non-strange light quarks: 

Ki(1270) = Kzp^ sin6' + Kip^ cos 9, 

i^i(1400) = Ksp^ cose- Kip^ sine, (2.26) 



with e PS -58° as implied from the study of D ^ i^i(1270)7r, i^i(1400)7r decays |39|. We use 
/ki(1270) = 175 MeV m to fix /3k{^p^) - Pk(^p,) = 0-303 GeV and obtain /^^^(uoo) = "87 MeV. 
Note that these /3j^(3p^), /3^(ip^) are close to 13k*- For the masses of Kip^ and Ksp^, we follow \24\ 
to determine them from the mass relations 2m\_^^ = "i^j(x232) ~'~"^hi{i380) "^^sp ~ "^A'i(i270) 
"^^"1(1400) ~ "^Jfi • -^^^ ^ ^'^'^ ^ systems, it is clear from Table Hill that |/^3/2| fs < fA^f^^ 
accordance with the expectation from HQS [cf. Eq. H2.7|) ]. 



III. COVARIANT MODEL ANALYSIS OF FORM FACTORS 

In this section we first review the analysis of the form factors for s-wave mesons within the 
framework of the covariant light-front quark model jlo| and then extend it to the p-wave meson 
case followed by numerical results and discussion. 

A. Form factors 



Form factors for P ^ P,V transitions are defined by 

{PiP")\V,\P{P')) = P.Uiq') + q,f-iq') 
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{V{P",e")\V^\P{P')) = e^^^pe'^-py^giq^), 

{V{P\e")\A,\P[P')) = -i{e";f(q^)+e*" ■P[p,a+[q^) + q,a^{q^)]}, (3.1) 

where P = P' + P" , q = P' — P" and the convention eoi23 = 1 is adopted. These form factors are 
related to the commonly used Bauer-Stech-Wirbel (BSW) form factors via 

„2 



Fne) = U{q% F--{q^) = U{q^) + -^f-{q% 



V^^'iq') = -{M' + M")g{q% (g') 



PV /2\ 



M' + M" ' 



A^^iq') = {M' + M") a^{q% 71^(^2) - 
where the latter form factors are defined by [4o| 



2M 



77 



{Y{P\^')\V^\P{P')) = 

{V{P",e")\A,\P{P')) = ,{(M' + M"K;Mr(g2)--£_ ^ 



- P A^^(a'^^ 



2M"'^q,[A7{e)-A^,''{q% 



with ^(0) = F(f^(0), AP^{<d) = vlf ^(0), and 



A^ (a ) 

The general expressions for P to low-lying p-wave meson transitions are given by |0] 



{Ay\P",e")\V, 
{A'/\P",e")\A, 
{A'/\P",e")\V, 
{A^/\P",e")\A, 
{T{P",e")\V, 
{T{P",e")\A, 



P{P')) = i[u+iq^)P^+u^{q^)q^], 

PiP')) = i{h/2{q'K+e"*-P[P,c'f{q') + q,c]/\q')]}, 

P{P')) = -qi,2{q^)e,,ape"*''P^qf', 

P{P')) = i[i,/2{q')e";+e"*-P[P,4\q') + q,c?!\q^)]}, 

P{P')) = -q3/2{q')e^ua(3e"*''P''qf', 

P{P')) = h{q^)e^,^f3e"*''^PxP''q^ 

P{P')) = -^{k{q')e%P^ + e'::pP'^pf^[P,b^{q')+q,b.{q')]}. 

^1/2(3/2) 1/2(3/2) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



The form factors ^1/2(3/2)1 c'^^^'^^^\ ^v^i^^/^J ^1/2(3/2) are defined for the transitions to the heavy 

1/2 3/2 

P^ {Pi ) state. For transitions to light axial- vector mesons, it is more appropriate to employ the 
L — S coupled states ^Pi and ^Pi denoted by the particles ^A and ^A in our notation. The relation 
between P^^'^^P'^^'^ and ^Pi, ^Pi states is given by Eq. (|2.6|) . The corresponding form factors 
and qiA{^A) foi' P ^ ^A i^A) transitions can be defined in an analogous way.^ 



The form factors i^A(^A)iC^^ '^\c^^ and qiA(^A) are dubbed as ^(w), c+(s+), c_ (s_) and q{r), respec- 
tively, in the ISGW model Q. 
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Note that only the form factors m+(5^), ii_ (g^) and ki^q^) in the above parametrization are 
dimensionless. It is thus convenient to define dimensionless form factors by^ 



{S{P")\A,\P{P')) 



' M'2-M"2 \^ps,2^ M'^-M"^ r^PS^ 2, 



{AiP",e")\V,\PiP')) = -iUmp - mA)e;vf%') - ' ^' P^V^'^q^) 

I mp — vnj^ 



e* -P' 
I2 



2mA-^q, Vi\q^) - V^\q^) 



{A{P\e")\A^\P(P')) = 6^,p.e*'^P^gM^^(52)^ (3_g) 



1 

mp — niA 

with 



and V/^(0) = Vo^^(O). They are related to the form factors in H3.3|) via 
Ff^iq') = -u^{q% FP{q^)=-u^{q^)-^u.{q% 
A^^iq') = -iM'-M")q{q'), V^iq') - 



M' - M" ' 



Vf^iq') = {W - M") c^{q% ¥1%') - V,^\q') = ^ c^{q\ (3.8) 

In above equations, the axial- vector meson A stands for A^/'^ or A^/'^. Besides the dimensionless 
form factors, this parametrization has the advantage that the q^ dependence of the form factors is 
governed by the resonances of the same spin, for instance, the momentum dependence of Fo^q^) is 
determined by scalar resonances. 

To obtain the P ^ M transition form factors with M being a ground-state s-wave meson or a 
low-lying p-wave meson, we shall consider the matrix elements 

{M{P")\V,-A,\P{P'))^BP^^, (3.9) 



where the corresponding Feynman diagram is shown in Fig. 1(b). We follow [10| to obtain P ^ P,V 
form factors before extending the formalism to the p-wave meson case. As we shall see, the P 
S {A) transition form factors can be easily obtained by some suitable modifications on P — > P (V) 
ones, and we need some extension of the analysis in [li| to the P ^ T case. 
For the case of M = P, it is straightforward to obtain 



„PP _ .3 f ,4 / H'pHp pp 



^ The definition here for dimensionless P ^ A transition form factors differs than Eq. (3.17) of where 
the coefficients (jnp ± m^) are replaced by (mp =F m^)- It will become clear in Sec. IV that this definition 
will lead to HQS relations for B Dq, Di transitions [cf. Eq. H4.7|l ] similar to that for B D, D* ones. 
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where 



qPP 



Tr[75(/i' + mlh^{ji[ + m'j75(- + m-2)], 



(3.11) 



Ni = p'l — nil + the subscript of S";/ stands for the transition vector current. As noted in 

the Introduction we consider the = frame As in the case, the p'{' integration picks 
up the residue p2 = P2 and leads to 



(") 



{")2n 



H 



/{//) 



M 

-II 

M 



N[N'^N2 



ttI ttII q 



PM 



'I'M ' 
II 



-it: 



— ilplljkjfD , 



where 



M, 



„2 _ + mf ^ pf + ml 



(3.12) 



(3.13) 



Xi X2 

with p'j_ = p'j_ — X2q±- In general, after the integration in B^^\ p'l can be expressed in terms 
of three external vectors, P' , q and a). Furthermore, the inclusion of the zero mode contribution 
cancels away the u dependence and in practice for p[ and N2 in 5^*^ under the integration, we 
have 3l 



»/ 

PlaPluPla 



N2 
Pl,N2 

P\pS>\vP\o,^2 



(3) 



_p p p /I (3) 



{Pf^Puqa + Pf^QuPa + q,.PuPa)A'i 
+ {q^iqi'Pa + qfiPuqa + P^,qyqa)Af^ + g^gi^^a^g^^ 
Z2, Xi7V2 ^ 0, 



(3) 



4^^Z2 + 



■P 



(2) 



'AfZ2 + 2^-^A^^Af^ 



+qt,quq 



3g2 



.(1) /.(3) 
^2 ^2 



1 

3^^ 



(2)n2 



(3.14) 



where A^j \ Z2 are functions of xi^2) p'_L) Pj_ ' q± and (;^, and their explicit expressions are given in 
[lfl| . We do not show the spurious contributions in the above equation since they vanish either after 
applying the above rules or after integration. The last rule on Pif^Pi^Pia^^ in the above equation, 
which is needed in the P ^ T calculation, is extended in this work. One needs to consider the 
product of four p'^'s. For completeness, the formulas for the product of four p'^s and the expressions 

(i) 

for A) ', Z2 can be found in Appendix B. 
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From Eqs. ^J^-^J^ one can obtain the form factors f±{q'^) for = -q^ < [see Eq. (|B3)) ]. 
We will return to the issue of the momentum dependence of form factors in the next sub-section. 
The explicit expressions for f± can be evaluated readily by using the explicit representations of 
N['-"\ given in Eqs. (1X1^ and dTTTl) . At = 0, the form factor /+(0) is reduced to the 
familiar form la. l4lll 



/+(0) 



1 



167r3 



dxd p'j_ ip"*{x,p'j_)ip'{x,p'^_) 



A' A" + p'l 



where 



/A'^+p'l^A"^+p'l 
A' = m'lX + 7712(1 — x), A" = mix + 7712(1 — x), 



with x = X2- 

For the P — > transition amplitude, we have 

N, 



B 



PS 



(2vr)^ 



(3.15) 



(3.16) 



(3.17) 



with 



SaS = Tr[(-i)(/i' + ?ni')7M75(/i + i>2 + m2)] 



(3.18) 



Thus, the P ^ S transition form factors are related to f± by 

u± = -f±{m'i ^ -ml h"p ^ h'^). 



(3.19) 



To be specific, we give the explicit forms of u±{q'^) obtained in the covariant light-front model 

Nr. 



u+{q ) 



u-{q') 



167r3 



h' h" 
dx.d^p'^J^^ 

^X2N{N'^ 



xi{M[ 



12 



M, 



Ii2\ 



X2q 



167r3 



I I I \ ii\2 , I I \2i ( II X ^2 

/"dx2dV±— x\X2M''^ +p1 + m'im2 + {m'l + 777-2) (x27?t,'i + 2:17772) 
J X2N[N{ y 



M 



112 



X2{q^ + q-P) 



{x2 - xi)M'^ + 2xiMlf - 2{m[ - m2)im[ - m'l) 



(3.20) 



It is ready to evaluate these form factors by using the explicit expressions of and h. Numerical 
study of these form factors will be given in the next sub-section. 

We next turn to the P ^ V,A transition form factors. For the P ^ V transition, we have 



.3 



N[N1N2 



(3.21) 



where 



(S 



V 



S 



Tr 



1 



7^ - ^{P'i - P2)v ) (// + "Ii)(7m - 77.75)(/i + m'^h^{- ^2 + ^2) 



V 



. (3.22) 
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As will be seen later, this expression of S is also useful for the P ^ T calculation, and hence 
its explicit representation is included in Appendix B. By the aid of Eos . (I3.12() and H3.14() . it is 
straightforward to obtain the P — > F form factors, g{q'^), f{q^)-, 0'±{<f ) 10]. For reader's conve- 
nience, the explicit forms of these form factors are summarized in Appendix B. Note that the vector 
form factor V{q^ = 0) = -{M' + M")g{q^ = 0) is consistent with that in 0, obtained by a 
Taylor expansion of the hy/N" term in g{q'^) [see Eq. ()B4|) ] with respect to p"^. To show this, we 
write 

d 

= — — 2 XoTl, -Oil- 



,^-2x2p'^-q±[- 



(3.23) 



and see that the second term on the right-hand side is needed when considering the — > limit 
of the p' ■ q±/q^ term in the integrand of g{q'^), while ©(x^Q'^) terms in the above equation vanish 
in the same limit. We perform the angular integration in the plane before taking the q^ ^ 



limit. After these steps, we obtain the same expression of V{q'^ = 0) as in |7I. I41||. 

The extension to P — > j4 transitions is straightforward and, as we shall see shortly, the resulting 
form factors have very similar expressions as that in the above case. For the P — > ^^4, transitions, 
we have 



B 



.3 



(2vr)^ 



d pi- 



N[N'{N2 
^1 N[N^N2 



qP^A ll*U 



qP ^A ll*V 
'-'fill ^ ) 



(3.24) 



where 



P3A 



P^A 



fcP^A 



gP^A- 



flV 



Tr 



Tr 



-(Pi 



■ P2)uj 75(/l + m'Dilfi - 7m75)(/i + "l'l)75(- h + "12 

71. - y^ip'i - P2)y ) (// - "^'i')(7m75 - 7m)(/i + K)75(- ^2 + m2 



^A 



( qP^A 



Tr 



qP^A- 
^A . 



Tr 



— {Pl - P2)u (// - "^i)(7m75 - 7/J(/l + "T''l)75(- h + ^2) 



— (pi - P2)y 75(/i + "ii)(7M - 7a.75)(/i + rn^)lb{- h + "^2) 



(3.25) 



We therefore have 5, 



^P'^A.P^A _ 

V(A) - ^A(y) 
only the XjW" terms are kept in S^^^. Consequently, 



S^y\ with the replacement —m", Wy 



W^AM- 



Note that 



q^'\q') 



f{q^) with {ml - 
g{q^) with (ml — 
aj^{(f') with {m'[ 



a-{q^) with 



[■m^ 



II ill 

-nil, '^v ~ 

m" h" 
-mi, Uy — 

-m'l, h'{, 

II 1,11 
-mi, hy 



uii 



'iAM 



II 

Wy 



h'l 



3A lA' 



h" 



ill II 

h. 



II \ 

■ W'^a^a) 

II \ 



3A lA' ""^3^ iaJ> 



(3.26) 
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where only the 1/W" terms in P — > form factors are kept. It should be cautious that the 
replacement of m'( — > — m" should not be applied to m" in w" and h" . These form factors can be 
expressed in the P^^'^ and P^^"^ basis by using Eq. ()2.6p . 
Finally we turn to the P ^ T transition given by 



PT 



.3 



(3.27) 



where 



cPT 



.^A-q+p[)x. (3.28) 

The contribution from the {—q)\ part is trivial, since q\ can be taken out from the integration, 
which is already done in the P ^ V case. Contributions from the Sj^^p'^^^ part can be worked out 
by using Eq. (|3.14|) . In particular, the calculation of k{q'^) and b-.{q^) needs to use the p'iPip'iN2 
formula. Putting all these together leads to 

Nr. 



b^iq' 



-9{q') 



+ 



167r3 



9/?' h" 

ax2d p_^- 



-X2N[N'{ 



(m;-m;')(4'^ + 4'0 



,(2)^ 



+ (^;' + _ 2m2)(Af + - m'Mi^ + A\'^) + -^{2Af^ + - A'i 



,(2) 



(1) 



,(3) 



,(3) 



+ 



K',^h'4. 167r3 



-m[{M"^ - N" - m 



h' h" 



X2N[N'{ 
"2 ml) - m'[(M'^ - N[ 



2(4'^+4'^)[m2(g'-iV(-iV(' 
'2 ml) - 2m[m';m2] 



mi 



Ii2\ 



mv 



+2(m; + m;')U2 ^2 + 



-4^ 

^2 "^1 



+ 16(m2 - m'i)(Ap^ + A^2^) + 4(2m; 
+4r ( [^^" + - 9' + 2(m; - m2)(m;' + m2)](2Af ^ + 2Af - A^^) 



mi 



m2)^l- 



(2) 



-a+{q^ 



+ 

V ' ""T 
(1) I /l(l) 



Nr. 



+ 2^f)Z2 



167r3 



h' h" 
r ,2 I ''■p'i'T 

ax2d p^ 



^4^|8(m2 - m;)(4^) + 2Af + 4^^) 



+- 



-2m;(4'^ + 4'0 + 4(2m; 
2 



m'/ 



m2)(4'^ + 4'^) + 2(m; + m'{){Al'^ + 2^'^ + A^ 



,(2) 



(2) 



,(2) 



1(2)^ 



2[M'2 + M"2 - g2 + 2(m; - m2)(m;' + m2)](4'^ + + ^5 

(1) 



.(3) 



(3) 



A 



(2) 



(2), 



+ [^2 _ ^/ _ _ ^ m;')'](^?^ + 24'^ + Af - A^l^ - A 



1(2) 



,(2) 



t(l) 



+ 



Nr 



h'^^^h'l 167r3 



/i' h" 
dxod^p'i — ^ ■ 

^X2iV(iV{' 



8(m2-m;)(4'^ + 24'^+4'^) 



(3) 



,(3)^ 



-6m; (.4f) + 4^^) + 4(2m; - m'/ - m2)(4^^ + 4^^) 
+2(3m; + m'l - 2m2)(^f ^ + 24^^ + 4^^) 



+- 



w. 



2[M'2 + M"2 - g2 + 2{m[ - m2)(m![ + m2)](4'^ + 2^?^ + 4'^ " 4'^ " ^f) 
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+2Z,(3Af - 24^) - + 2^ (64') - 64^)4'^ + ^{A?f - Af^) 
+ - 2M'2 + N[- N'{ - (m; + m'if + 2(m; - ma)^] 



X (4^) + 24^) +42) -4^) -4^)^ 



(3.29) 



To summarize, equipped with the exphcit expressions of the form factors /+(g^), /_(g^) 
[Eq. ((HSl)] for P ^ P transitions, u+{q^),U-{q^) [Eq. (jS^Dl)] for P ^ 5 transitions, 
9{q^)J{q^),a+{q^),a^{q^) [Eq. (jEl)] for P ^ F transitions, £(g2), g(g2)^ c+(g2), c_ (^2) [^^^ ^^SHl)] 

for P ^ A transition and hi^q"^), k(q^),b+{q'^),b^{q'^) [Eq. 1)3.291) ] for P ^ T transitions, we are 
ready to perform numerical studies of them. The P ^ S, A,T transition form factors are the main 
new results in this work. 



B. Comments on zero-mode effects 

In the present paper we have followed and extended the work of Jaus to the p-wave meson 
case. As stressed by Jaus, there are two classes of form factors and decay constants. There is one 
class of form factors like Fi{q^) for transitions between pseudoscalar mesons, V{(f') and A2{q^) 
for transitions between pseudoscalar and vector mesons, and the pseudoscalar decay constant fp 
that are free of zero mode contributions. Another class of form factors like Ai{q^) (or f{q^)) and 
the vector decay constant fv are associated with zero modes. The full vector vertex operator for 
^5i-state meson has the expression (see Table I) 

iHv[l^.-^{Pl-P2)^^. (3.30) 

To begin with, we first consider the "simple" vector meson vertex without the 1/W part in the 
above expression. Jaus employed a simple multipole ansatz for the meson vertex function 

as the starting point of his simple covariant toy model. Then the zero mode contributions can be 
systematically calculated in this toy model. Note that the vertex function 1)3. 31() is not symmetric 
in the four momenta of the constituent quarks and hence can hardly be considered a realistic 
approximation of the meson vertex. To remedy this difficulty, Bakker, Choi and Ji (BCJ) 0] 
proposed to replace the point gauge-boson vertex 7^(1 — 75) by 

It is easily seen that the two methods due to Jaus and BCJ should give the same result for form 
factors, but may lead to different results for decay constants. Indeed, Eq. (3.9) (without the 1/W 
part) of Jaus 'l4j for the form factor f{q^) agrees with Eqs. (37) and (38) of BCJ 11 j. Moreover, 



it is interesting to notice that Eq. (3.16) of Jaus Q] for the decay constant fv (by considering the 
longitudinal polarization case as in BCJ) also agrees with Eq. (41) of BCJ though the analytic 
expressions for the respective vertex functions Hy are different [4j|. Since the associated trace for 
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fv is free of minus components of the internal momenta and there are no zero modes in that case. 
Therefore, to the level without the 1/W part in the vertex operator (|3.3U() . there is no discrepancy 
between Jaus and BCJ and both fv as well as f{q'^) are free of zero mode effects. 

However, things are very different when the full vector meson vertex (|3.3Uj) is used. The 1/W 
part of the trace contains minus components of the momenta and the zero mode problem must 
be faced. While BCJ claimed that both fv and /(g^) are immune to the zero mode even for 
the full vector meson vertex, Jaus obtained non-trivial zero mode contributions. It appears to us 
that the controversy about the role played by the zero mode lies in the fact that Jaus and BCJ 
have different procedures for identifying zero-mode effects. In the covariant light-front approach 
of Jaus the decomposition of the current-induced matrix element into 4- vectors will require 
to introduce a lightlike 4-vector a; which is not covariant. Zero modes are required to eliminate 
the spurious Co dependence. BCJ decompose the propagator into on-shell and instantaneous (not 
on-shell) parts and show that only the latter part can be the origin of a zero mode contribution. 
More precisely, the contour integration over p~ in Jaus is not a regularized one, while in BCJ the 
contour can be regularly closed due to the presence of the non-local boson vertex and the zero 
modes display their effects at the level of p'^ (see Sec. II. B. 2 of 0] for more detail about the BCJ 
approach for zero modes). 

The covariant toy model cannot be generalized beyond the simple meson vertex given by 
namely, there are some possible residual Co contributions. In Jaus has developed a method 
that permits the calculation of the contribution of zero modes associated with the current-induced 
matrix element. Through the study of the angular condition imposed on helicity amplitudes, several 
consistency conditions can be derived under some plausible assumptions and used to determine the 
zero mode contributions. Within this approach, both fv and /(g^) receive additional residual 
contributions (see Eqs. (3.16) and (3.9) of respectively) which can be expressed in terms of 
-BI"^^ and Cn™^ functions defined in Appendix B.^ These functions depend on p'{' and behavior 
like {p']~y{p'i~y ■ Jaus then gave a counting rule for detecting zero modes For the B functions 
i < j, there is no zero mode contribution and the value of -Si"^^ can be calculated unambiguously 
at the spectator quark pole. For the C functions i > j + 1 and the value of Cn™'* is the sum of 
a spectator quark pole term and an unknown zero mode contribution. These B and C functions 
vanish in the covariant toy model, as it should be. Beyond the toy model, the C functions contain 
unknown zero-mode contributions. Jaus used some consistency conditions to fix some of the C 
functions. 

Several remarks are in order, (i) Our meson light-front vertex functions 1)2. are symmetric 
in quark momenta. However, B and C functions do not appear in fv [Eq. 1)2.22(1 ] and f{q'^) [Eq. 
()B4|) ] for two reasons. First of all, we have ^ = from Eqs. (Hnil) and (jUl). Second, we contract 
Aj^ [see Eq. (|2.2U() ] and Sj^^ [Eq. (|3.22|) ] with the transverse polarization vector e^(ib). We have 
checked explicitly that for the vertex functions given in (|2.11() , the coefficients i? j*^ under integration 



^ The remaining spurious uj contribution to the form factor a^(q'^) cannot be determined in the same 
manner. 
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FIG. 2: Form factors F\{c^) and Fo((?^) for i? ^ D and B ^ Dq transitions. 



(see Appendix B) are numerically almost vanishing and the form factor f{q^) is affected at most at 
one percent level. For our purposes, we can therefore neglect all the residual contributions to the 
form factors, (ii) The derivation of the decay constant fp and the form factors f+{q'^),g{q'^),a+{q'^), 
u+{q'^),q{q'^),c^{q'^), and h{q'^),b^{q'^) does not depend on iV2 and those relations connected to N2 
(see Eq. (j3.14|) and recall that N2 = Z2 — c[^^). These form factors are free of zero mode effects 
and can be obtained using the conventional light-front approach, (iii) Zero mode effects vanish in 
the heavy quark limit (see Sec. IV). For example, the HQS relation fp = fy indicates that fv is 
immune to the zero mode contribution. 



C. Form- factor momentum dependence and numerical results 

Because of the condition q~^ = we have imposed during the course of calculation, form factors 
are known only for spacelike momentum transfer q^ = —q'j_ < 0, whereas only the timelike form 
factors are relevant for the physical decay processes. It has been proposed in 6] to recast the form 
factors as explicit functions of in the spacelike region and then analytically continue them to the 
timelike region. Another approach is to construct a double spectral representation for form factors 



at < and then analytically continue it to > region 4^. It has been shown recently that, 
within a specific model, form factors obtained directly from the timelike region (with 0+ > 0) are 
identical to the ones obtained by the analytic continuation from the spacelike region [ll|. 

In principle, form factors at > can be evaluated directly in the frame where the momentum 
transfer is purely longitudinal, i.e., q± = 0, so that q^ = q~^q~ covers the entire range of momentum 
transfer 0]. The price one has to pay is that, besides the conventional valence-quark contribution, 
one must also consider the non- valence configuration (or the so-called Z-graph) arising from quark- 
pair creation from the vacuum. However, a reliable way of estimating the Z-graph contribution is 



still lacking unless one works in a specific model, for example, the one advocated in Fortunately, 
this additional non-valence contribution vanishes in the frame where the momentum transfer is 
purely transverse i.e., q^ = 0. 

To proceed we find that except for the form factor V2 to be discussed below, the momentum 
dependence of form factors in the spacelike region can be well parameterized and reproduced in the 
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FIG. 3: Form factors V{ct), ^o(g^), Mist) and ^2(0'^) for B ^ D* transitions. 



three-parameter form: 



m 



(3.33) 



^B{D)>+K(fl'T^B(D)y ' 

for B{D) M transitions. The parameters o, b and F{0) are first determined in the spacehke 
region. We then employ this parametrization to determine the physical form factors at > 0. In 
practice, the parameters a, b and -F(O) are obtained by performing a 3-parameter fit to the form 
factors in the range -20 GeV^ < g2 < for 5 decays and -10 GeV^ < q"^ < iov D decays. These 
parameters are generally insensitive to the range to be fitted except for the form factor V2(g^) in 
B(D) ^Pi,P^^'^ transitions. The corresponding parameters a and b are rather sensitive to the 
chosen range for g^. This sensitivity is attributed to the fact that the form factor V2((?^) approaches 



TABLE IV: Form factors oi D ^ tt, p, oq (1450), oi (1260), 61 (1235), 02 (1320) transitions obtained 
in the covariant light-front model are fitted to the 3-parameter form Eq. (|3.33|) except for the 
form factor V2 denoted by * for which the fit formula Eq. (|3.34|) is used. All the form factors are 
dimensionless except for h, 6+, b^ with dimensions GeV~^. For the parameter Pt appearing in the 
tensor-meson wave function, we assume that it is the same as the P parameter of p-wave meson 
with the same quark content. 



F 


m 


^(^max) 


a 


b 


F 


m 


-^(^max) 


a 


b 


pDn 


0.67 


2.71 


1.19 


0.36 


^0 


0.67 


1.16 


0.50 


0.01 


yDp 


0.86 


1.36 


1.24 


0.48 




0.64 


0.93 


1.07 


0.54 




0.58 


0.71 


0.51 


0.03 




0.48 


0.68 


0.95 


0.30 


pDao 


0.52 


0.54 


1.07 


0.26 


rpDao 
^0 


0.52 


0.52 


-0.08 


0.03 




0.20 


0.22 


0.98 


0.20 




0.31 


0.34 


0.85 


0.49 


yDai 


1.54 


1.53 


-0.05 


0.05 




0.06 


0.06 


0.12 


0.10 




0.11 


0.13 


1.08 


0.54 




0.49 


0.54 


0.89 


0.28 




1.37 


1.45 


0.46 


0.05 


yDb, 


-0.10* 


-0.11* 


0.21* 


0.67* 


h 


0.188 


0.208 


1.21 


1.09 


k 


0.340 


0.338 


-0.07 


0.12 




-0.084 


-0.091 


0.97 


0.58 


b- 


0.120 


0.133 


1.15 


0.66 
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FIG. 4: Form factors A{q^), Vb(g2), and V2(g^) for B D{' ' ' transitions. 
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to zero at very large — \q^\ where the three-parameter parametrization H3.33() becomes questionable. 
To overcome this difficulty, we will fit this form factor to the form 

F(0) 



(3.34) 



BKi, 



and achieve a substantial improvement. For example, we have a = 2.18 and h = 6.08 when 



TABLE V: Same as Table HVI except for D K, K*, K(;(1430), Kip^ , Ksp^ , K|(1430) transitions. 



F 


m 


-^('i'max) 


a 


b 


F 


m 


F{Qmax) 


a 


b 


pDK 


0.78 


1.57 


1.05 


0.23 




0.78 


0.99 


0.38 


0.00 


yDK* 


0.94 


1.33 


1.17 


0.42 




0.69 


0.92 


1.04 


0.44 




0.65 


0.75 


0.50 


0.02 


aDK* 
^2 


0.57 


0.75 


0.94 


0.27 




0.48 


0.51 


1.01 


0.24 


DK* 


0.48 


0.50 


-0.11 


0.02 




0.10 


0.11 


1.03 


0.48 




0.44 


0.47 


0.80 


0.27 




1.53 


1.58 


0.39 


0.05 


DKip 
^2 ' 


-0.09* 


-0.09* 


-0.16* 


0.51* 




0.98 


1.05 


0.92 


0.17 




0.34 


0.38 


1.44 


0.15 




2.02 


2.02 


-0.01 


0.03 




0.03 


0.03 


-0.18 


0.10 


h 


0.192 


0.205 


1.17 


0.99 


k 


0.368 


0.367 


-0.04 


0.11 




-0.096 


-0.102 


1.05 


0.58 


6_ 


0.137 


0.147 


1.17 


0.69 
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FIG. 5: Form factors k{q^), h{q^), h+{q^) and b_{q^) for B ^ D2 transitions. Except for the 
dimensionless k{q'^), all other form factors are in units of GeV~^. 



is fitted to Eq. (|3.3ci|) and they become a = 1.78 and 6 = 2.12 (see Table rVTI]) when the fit formula 
Eq. H3.34|) is employed. 

In Tables ITVl - IVIIII we show the form factors and their dependence for the transitions B{D) 
7r,/),ao(1450),ai (1260), 61 (1235), 02 (1320), B{D) K, K*, K(;( 1430), Kip^, Ksp^, 1430) and 
B D,D*,L>i5 (2308), ^)l/^L>f^L>| (2460). The 6 ^ c transition form factors are plotted in 
Figs. Because the quark contents of the /o, /i, /2 mesons lying in the mass region of 1.3 — 1.7 

GeV are not well established, we will not consider them in this work. In calculations, we have taken 
.he meson masses from E- The masses of D'„ and D, have been measured recently by Belle to 
be 2308 ± 17 ± 15 ± 28 MeV and 2427 ± 26 ± 20 ± 15 MeV, respectively Q. Since L>i(2427) and 
Z)'^(2420) are almost degenerate, we shall take m 1/2 ~ m 3/2 ~ 2427 MeV. 

Several remarks are in order: 

1. Many form factors contain terms like {p'±-q±)/q'^ in their integrands. At first sight, it appears 



TABLE VI: Same as Table HVl except for B vr, p, ao (1450), d (1260), 61 (1235), 02 (1320) transi- 



tions. 



F 


m 


i^(dax) 


a 


b 


F 


m 


-^(^max) 


a 


b 




0.25 


1.16 


1.73 


0.95 


TpBir 


0.25 


0.86 


0.84 


0.10 




0.27 


0.79 


1.84 


1.28 




0.28 


0.76 


1.73 


1.20 




0.22 


0.53 


0.95 


0.21 




0.20 


0.57 


1.65 


1.05 


pBao 


0.26 


0.68 


1.57 


0.70 


T-i_Bao 
-^0 


0.26 


0.35 


0.55 


0.03 


j^Bai 


0.25 


0.76 


1.51 


0.64 


^0 


0.13 


0.32 


1.71 


1.23 


yBai 


0.37 


0.42 


0.29 


0.14 




0.18 


0.36 


1.14 


0.49 


^Bb^ 


0.10 


0.23 


1.92 


1.62 




0.39 


0.98 


1.41 


0.66 


yBb, 


0.18 


0.36 


1.03 


0.32 




-0.03* 


-0.15* 


2.13* 


2.39* 


h 


0.008 


0.015 


2.20 


2.30 


k 


0.031 


0.010 


-2.47 


2.47 




-0.005 


-0.011 


1.95 


1.80 




0.0016 


0.0011 


-0.23 


1.18 
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TABLE VII: Same as Table |IV| except for B 




, K'ip^ , K2 (1430) transitions. 


F 


m 




a 


b 


F 


F(0) 


-^(^max) 


a 


b 


pBK 


0.35 


2.17 


1.58 


0.68 


Fi^ 


0.35 


0.80 


0.71 


0.04 


yBK 


0.31 


0.96 


1.79 


1.18 




0.31 


0.87 


1.68 


1.08 


aBK* 
^1 




58 


93 


1 9 


^2 


0.24 


70 




98 






n 70 

U. / u 




U.U'4: 


^0 




u.oo 




U.U(J 








1 47 






1 4 




1.62 


1.14 




0.39 


0.42 


0.21 


0.16 




0.17 


0.30 


1.02 


0.45 




0.11 


0.25 


1.88 


1.53 


BKip 


0.41 


0.99 


1.40 


0.64 




0.19 


0.35 


0.96 


0.30 




-0.05* 


-0.16* 


1.78* 


2.12* 


h 


0.008 


0.018 


2.17 


2.22 


k 


0.015 


0.004 


-3.70 


1.78 




-0.006 


-0.013 


1.96 


1.79 




0.002 


0.002 


0.38 


0.92 


TABLE VIII: Same 


as Table IIVI except for B 


^ D, D* 


did\'\d 


\^'^,D2 transitions. 


For the 


purpose 


of comparin 


g with heavy quark symmetry, the form factors u±,c±,i,q are 


also shown. 


F 


m 


-^(^max) 


a 


b 


F 


m 


-^(^max) 


a 


b 


pBD 


0.67 


1.22 


1.25 


0.39 


TpBD 


0.67 


0.92 


0.65 


0.00 


yBD' 


0.75 


1.32 


1.29 


0.45 


aBD* 
^0 


0.64 


1.17 


1.30 


0.31 


ABD' 


0.63 


0.83 


0.65 


0.02 


aBD* 
^2 


0.61 


0.95 


1.14 


0.52 


pBDt, 


0.24 


0.34 


1.03 


0.27 




0.24 


0.20 


-0.49 


0.35 




-0.12 


-0.14 


0.71 


0.18 


BD^'^ 


0.08 


0.13 


1.28 


-0.29 


BD^^^ 

vf ' 


-0.19 


-0.13 


-1.25 


0.97 


T/ 1 


-0.12 


-0.14 


0.67 


0.20 


ABD{'^ 


0.23 


0.33 


1.17 


0.39 


bdI^^ 


0.47 


0.70 


1.17 


0.03 


^bdT 


0.55 


0.51 


-0.19 


0.27 


bd\'^ 

^2 


-0.09* 


-0.17* 


2.14* 


4.21* 




-0.24 


-0.34 


1.03 


0.27 


n_ 


0.31 


0.42 


0.86 


0.20 


h/2 


0.56 


0.38 


-1.25 


0.97 


91/2 


0.041 


0.050 


0.71 


0.18 


1/2 

c+ 


-0.042 


-0.050 


0.67 


0.20 




0.045 


0.055 


0.71 


0.20 


4/2 


-1.56 


-1.45 


-0.19 


0.27 


93/2 


-0.079 


-0.114 


1.17 


0.39 


cf 


-0.032* 


-0.061* 


2.14* 


4.21* 


3/2 
C_ 


-0.027 


-0.026 


0.03 


0.45 


h 


0.015 


0.024 


1.67 


1.20 


k 


0.79 


1.12 


1.29 


0.93 




-0.013 


-0.021 


1.68 


0.98 


6_ 


0.011 


0.016 


1.50 


0.91 



that linear p'^ terms will not make contributions after integrating over p'^. But this is not 
the case. As noted before, a Taylor expansion of the h"j^,.jNi term with respect to p"^ will 
generate a term proportional to • (7_L [cf. Eq. ()3.23|) ]. When combined with the •q'j,)/^^ 
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term in the integrand of transition form factors, this leads to 



d Pi_ 2 = -n dp_^pj_ (3.35) 



2 

in the = frame. In analytic studies, however, it is more convenient to utilize the identity 
obtained in 0] 

/ dx2 ciV±^^4^25f ) = / dx2 dV±-^4^ (x,Z2 - 2Af^) = 0. (3.36) 

Using the expressions of Z2 and A^'^ given in Eq. (|B9|) , it is easily seen that the {p'j_ ■ qi_)/q^ 
term under integration can be related to other g-independent quantities. The above identity 
allows us to integrate out the p'^ ■ q± term without explicitly performing the Taylor expansion 
of h'^f/N^. Instead of using Eq. (|3.35() or (|3.36|) we have taken into account such effects in 
numerical calculations by substituting the relation p'j_ = p'j_ — X2q± into h'j^j/Ni . 

2. Owing to the less energy release, form factors for D ^ tt, p, ■ ■ ■ and D ^ K, K* , K** transi- 
tions are more sensitive to the masses of charmed and light quarks. For this we can utilize the 
form-factor ratios ry = V^"^ {0)/Af^ (0) and r2 = A^^ (0) / Af^ {0) measured in D ^ V£u 
decays to constrain the quark masses. The most recent and most precise measurement of 
D+ K*^i+v by FOCUS yields H 

rv{D K*) = 1.504 ± 0.057 ± 0.039 , r2{D K*) = 0.875 ± 0.049 ± 0.064 . (3.37) 

The best quark masses mu,ms and rric obtained in this manner are listed in Eq. (|2.24jl . 
Using this set of quark masses and fixing f3^ = 0.3073 GeV from = 237 MeV we have also 
computed the similar form factor ratios for (f) transition and found 

rv{Ds -^(j)) = 1.569 , r2{Ds (p) = 0.865 , (3.38) 

in goo d ag reement with the very recent FOCUS measurement of the Df — > (pfi^v form factor 
ratios |45| 

rv/(-Ds ^ 0) = 1-549 ±0.250 ±0.145, r2(-Ds ^ </>) = 0.713 ± 0.202 ± 0.266 . (3.39) 

3. In the absence of any information for the parameter (3t appearing in the wave function of 
tensor mesons, we have taken /3t to be the same as the (5 parameter of the p-wave meson 
with the same quark content, for example, f3{D2) = P{Dq) = 0.331. Note that among the 
four P ^ T transition form factors, the one k{q'^) is particularly sensitive to fSx- It is not 
clear to us if the complicated analytic expression for k{q'^) in Eq. (|3.29|) is not complete. To 
overcome this difficulty, we apply the heavy quark symmetry relations in Eq. (|4.4|) below to 
obtain k{q'^) for B D2 transition 



k{q ) = mBmn; 1 + 



hiq')-lb+{q') + lb^{q') 



2mBmD* 

This can be tested in D^tt" decays to be discussed below in Sec. III.E. 



(3.40) 
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4. For heavy-to-heavy transitions such as i? — > D,D*,D**, the sign of various form fac- 
tors can be checked by heavy quark symmetry. In the heavy quark hmit, heavy quark 

1/2 

symmetry requires that the form factors n_, ^i/2) 9i/2i c_ ,h,k and 6_ be positive, while 
.£3/2, ^3/2) c+ ,c/ , c_ and 6+ be negative [see Eqs. H4.3() — ()4.5() ]. Our results are indeed 
in accordance with HQS (see Table IVlTTj) . 

5. For P ^ A transitions, the form factor Vq is always positive, while the sign of other form 
factors A, Vi, V2 depends on the process under consideration, for example, they are all positive 
in B{D) ai, Ksp^ transitions and negative in S — > Di transitions. 

6. The form factors of B to light axial- vector meson transitions obey the relations Vq""^ 



< 



A^^ < V,^"^ and V,^""''^ < A^* < V^^'^ . 



7. It is pointed out in 7] that for B ^ D, D* transitions, the form factors Fi, Aq, A2,V exhibit 
a dipole behavior, while Fq and Ai show a monopole dependence. According to the three- 
parameter parametrization 1)3. 33() . the dipole behavior corresponds to 6 = (a/2)^, while 6 = 
and a 7^ induces a monopole dependence. An inspection of Tables HVT - IVIIII indicates that 
form factors F^^, Af^*, Fq^, Fq""*^ and fI^^^'^^ have a monopole behavior, while 
yBD^^ jsBdT^ (^)^, ^^(^)^3p,^ y^^^..^ pBiD)ao pDK*, ^^^^ ^ ^^^^^ dependence. 

8. In the heavy quark limit, Ff^°{q^) = — q^/{mB - "Id^)^], while Ff^^q"^) = 

- q^/{mB + rno)'^] (see Eqs. (jOI) and This explains why Fi and Fq in 

the B Dq transition deviate at large faster than that in the B ^ D case (Fig. [2). 

9. Unlike the form factor Fq in P — > gqjKq transitions which is almost flat in its q'^ behavior, 
p^^o decreasing with q"^ as it must approach to zero at the maximum q"^ when mg 00 
[see Eq. ()4.7() ]. In general, form factors for P ^ S transitions increase slowly with 
compared to that for P ^ P ones. For example, F^"o(0) ~ F^'^(O) at q^ = 0, while at 
zero recoil ^^'"0(9^^^) < F^'^ (ql^.^^). Note that the form factors of B ^ qq 01 B ^ Kq 
are similar to that of 5 ^ vr or 5 ^ at = 0, while Fff^iO) <C Fff{0). This can be 
understood from the fact that P ^ S form factors are the same as P — > P ones except for the 
replacement of m'( — m" and hp —h'g [see Eq. (|3.20p ]. Consequently, the A!' term in 
Eq. I|3.16jl is subject to more suppression in heavy-to-heavy transitions than in heavy-to-light 
ones. We shall see in Sec. III.E that the suppression of the B Dq form factor relative to 
the B ^ D one is supported by experiment. 



10. To determine the physical form factors for B{D) E:i(1270), i^i (1400), B 

(2427), Di (2420), Dsi (2460), 1)51(2536) transitions, one needs to know the mixing angles 
of ^Pi -3 Pi [see Eq. (IT^ J and dI^"^ - D^^'^ . As noted in passing, the mixing angle for Ki 
systems is about —58° as implied from the study of D ^ -ftri(1270)7r. i^i(1400)7r decays js^. 
A mixing angle 9di = (5.7 it 2.4)° is obtained by Belle through a detailed B — > D*Tnr anal- 
ysis while 0Dsi ~ 7° is determined from the quark potential model 3^ as the present 
upper limits on the widths of Dsi{2A60) and 1)^1(2536) do not provide any constraints on 
the dIi^ — D^i^ mixing angle. 
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TABLE IX: Form factors of D ^ vr, p, K, K* transitions at (7^ = in various models. 



Model 


F/j7(0) .4^^(0) ^f^(O) ^^''(O) V^P{^) 


F/^o^(0) Ao^^*(0) A 


f^*(o) ^?^*(o) y^^*(o) 


This work 


0.67 


0.64 


0.58 


0.48 


0.86 


0.78 


0.69 


0.65 


0.57 0.94 


MS [46] 


0.69 


0.66 


0.59 


0.49 


0.90 


0.78 


0.76 


0.66 


0.49 1.03 


QSR [42 


0.5 


0.6 


0.5 


0.4 


1.0 


0.6 


0.4 


0.5 


0.6 1.1 


BSW [40] 


0.69 


0.67 


0.78 


0.92 


1.23 


0.76 


0.73 


0.88 


1.15 1.23 



D. Comparison with other model calculations 



It is useful to compare our results based on the covariant light-front model with other theoretical 
calculations. Except for the Isgur-Scora-Grinstein-Wise (ISGW) quark model all the existing 
studies on mesonic form factors focus mainly on the ground-state s-wave to s-wave transitions. 
For P ^ P,V form factors we choose the BSW model |4o| . the Melikhov-Stech (MS) model 
QCD sum rule (QSR) |47| and light-cone sum rules (LCSR) |43] for comparison. Shown in Tables 
IIXHXII are {D,B) n, p, K, K* , D, D* transition form factors calculated in various models. We 
see that the covariant light-front model predictions are most close to that of the MS model except 
for B — > K* transitions. 



ISGW model : 

Before our work, the ISGW quark model is the only model that can provide a systematical 
estimate of the transition of a ground-state s-wave meson to a low-lying p-wave meson. However, 
this model is based on the non-relativistic constituent quark picture. In general, the form factors 
evaluated in the original version of the ISGW model are reliable only at = g^, the maximum 
momentum transfer. The reason is that the form-factor q'^ dependence in the ISGW model is 
proportional to exp[— (g,^ — g^)] and hence the form factor decreases exponentially as a function of 
{Qm ~ ^^)- This has been improved in the ISGW2 model 30] in which the form factor has a more 
realistic behavior at large (g^ — (f') which is expressed in terms of a certain polynomial term. In 
addition to the form-factor momentum dependence, the ISGW2 model incorporates a number of 
improvements, such as the constraints imposed by heavy quark symmetry, hyper fine distortions of 
wave functions, • • •, etc. jiol. 

The ISGW2 model predictions for B D** transition form factors are shown in Table IXIII 



TABLE X: Form factors of B ^ tt, p, K, K* transitions at = in various models. 



Model 


^1^0" (0) ^^""(O) AfP{fd) A^P{{)) V^P{{)) 


F^^i^{Q) A^^\<d) Af^*(0) ^f^*(0) y^^*(0) 


This work 

MS [46] 
LCSR [48] 
BSW [40] 


0.25 0.28 0.22 0.20 0.27 
0.29 0.29 0.26 0.24 0.31 
0.26 0.37 0.26 0.22 0.34 
0.33 0.28 0.28 0.28 0.33 


0.35 0.31 0.26 0.24 0.31 
0.36 0.45 0.36 0.32 0.44 
0.34 0.47 0.34 0.28 0.46 
0.38 0.32 0.33 0.33 0.37 
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TABLE XI: Form factors oi B ^ D, D* transitions at g'^ = in various models. 



Model 




Ao^^*(0) 


ylf^*(0) 


ylf^*(0) 


F^^'(O) 


This work 


0.67 


0.64 


0.63 


0.62 


0.75 


MS [46] 


0.67 


0.69 


0.66 


0.62 


0.76 


BSW [40] 


0.69 


0.62 


0.65 


0.69 


0.71 



Note that form factors Fq °{q^), Vi ^ {q^) (or ^1/2(9^)) cannot be parameterized in the form of 

(|3.33jl or (|3.34|) since they vanish at certain g^, e.g. around « 8GeV^ for ^ iq^)- We see 
from the comparison of Table IXIII with Table IVIIII that (i) the form factors at small q^ obtained in 
the covariant light-front and ISGW2 models agree within 40%, and (ii) as q^ increases, Ff^°{q'^), 

3/2 BD'^^^ 

j4^^i (9^), ^ (9^)) h[q^), \bj^{q^)\ and h-[(f') increase more rapidly in the LF model than 

BD* BD^^^ 

those in the ISGW2 model, whereas Fq " (g-^) and \ V-^ ^ {q'^)\ decrease more sharply in the latter 
model so that they even flip a sign near the zero recoil point. 

The fact that both LF and ISGW2 models have similar B D** form factors at small q^ implies 
that relativistic effects could be mild in i? — > D** transitions. Nevertheless, relativistic effects may 
manifest in heavy-to-light transitions, especially at the maximum recoil. An example is provided 
shortly below. 

Others : 

Based on the light-cone sum rules, Chernyak j3jQ] has estimated the B ao(1450) transition 



TABLE XII: Form factors ol B ^ D** transitions calculated in the ISGW2 model. 



F 


m 




a 


b 


F 


m 


-^('i'max) 


a 


b 


pm 


0.18 


0.24 


0.28 


0.25 




0.18 


-0.008 






j^BdT 


-0.16 


-0.21 


0.87 


0.24 




0.18 


0.23 


0.89 


0.25 


^BD\^^ 


-0.19 


0.006 






bdY^ 

^2 


-0.18 


-0.24 


0.87 


0.24 




0.16 


0.19 


0.46 


0.065 


BD^/^ 
^0 


0.43 


0.51 


0.54 


0.074 


vf ^ 


0.40 


0.32 


-0.60 


1.15 


V2 ' 


-0.12 


-0.19 


1.45 


0.83 




-0.18 


-0.24 


0.88 


0.25 


U- 


0.46 


0.62 


0.87 


0.25 




0.54 


-0.016 






qi/2 


0.057 


0.074 


0.87 


0.24 


1/2 

c+ 


-0.064 


-0.083 


0.87 


0.24 


1/2 

cj 


0.068 


0.088 


0.87 


0.24 


h/2 


-1.15 


-0.90 


-0.60 


1.15 


93/2 


-0.057 


-0.066 


0.46 


0.065 


cf 


-0.043 


-0.066 


1.45 


0.83 


3/2 

C_ 


-0.018 


-0.013 


0.23 


5.38 


h 


0.011 


0.014 


0.86 


0.23 


k 


0.60 


0.68 


0.40 


0.68 


K 


-0.010 


-0.013 


0.86 


0.23 


6_ 


0.010 


0.013 


0.86 


0.23 



29 



form factor and obtained F^q''(0) = 0.46, while our result is 0.26 and is similar to the B ^ n 
form factor at = 0. For B ai(1260) form factors, there are two existing calcu ations: one in a 
,uark-me»„ .ode, (CQM) and the Che. based on the QCD sum rule (QSR) 0- 

The results 

are quite different, for example, Vq^"^(0) obtained in the quark-meson model, 1.20, is larger than 
the sum-rule prediction, —0.23 it 0.05 , by a factor of five apart from a sign difference. Predictions 
in various models are summarized in Table IXlTll If ai(1260) behaves as the scalar partner of the p 
meson, it is expected that Vf^""^ is similar to j4^^. Therefore, it appears to us that a magnitude of 
order unity for Vg^'*^ (0) as predicted by the ISGW2 model and CQM is very unlikely. Notice that 
the sign of the form factors predicted by QSR is opposite to ours. In hadronic B — > aiP decays, 
the relevant form factors are V^"^"^ and F^^^ under the factorization approximation. Presumably, 
the measurement of B^ a^vr" will enable us to test Vf^""^ . 



E. Comparison with experiment 



There are several experimentally measured decay modes, namely, B — > DD** and B^ D**Tr~ 
decays, which allow to test our model calculations of decay constants and form factors for p-wave 
charmed mesons D** and D** . 

For B DD** decays, they proceed only via external VF-emission and hence can be used to 
determine the decay constant of D** . More precisely, their factorizable amplitudes are simply given 
by 

A(B ^ DD*;) = ^V,tV:,ar{D*;\{sc)\0){D\{cbm, (3.41) 

where_{Qig2) = ft7/x(l — 75)92 and oi is a parameter of order unity. The recent Belle measurements 
read |37| 

B[B DD*o(2317)]i3[L'*o(2317) ^ L»,7r°] = {8.5tji ± 2.6) x 10"^, 

B[B DDsii2m)]B[Dsii2m) Z)*7r°] = (17.8t|^ ± 5.3) x 10"^. (3.42) 

The D*o(2317) width is dominated by its hadronic decay to Dgir^ as the upper limit on the 
ratio T{D*Q D*-f)/T{D*Q D^tt^) was set to be 0.059 recently by CLEO 16]. Therefore, 



TABLE XIII: B ai(1260) transition form factors at = in various models. The results of 
CQM and QSR have been rescaled according to the form- factor definition in Eq. 1)3. 6() . 



Model 


^Ba, (0) 


yo^'^^(o) 






This work 


0.25 


0.13 


0.37 


0.18 


ISGW2 [211] 


0.21 


1.01 


0.54 


-0.05 


CQM [49] 


0.09 


1.20 


1.32 


0.34 


QSR [50] 


-0.41 ±0.06 


-0.23 ±0.05 


-0.68 ±0.08 


-0.33 ± 0.03 
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0.94 < B[D*q{2317) Dsir"] < 1.0. It follows from Eqs. (jO^ and f^Tm that (see ^ for detail) 



/d:^ « 47 - 73 MeV, (3.43) 

for ai = 1.07. To estimate the branching ratios of D*-k^ in the -D5i(2460) decay, we need some 
experimental and theoretical inputs. There are two different measurements of the radiative mode 
by Belle: a value of 0.38 ± 0.11 ± 0.04 for the ratio Ds^/D*tt^ is determined from B DDgi 
decays while the result of 0.55 it 0.13 it 0.08 is obtained from the charm fragmentation of 
e~^e~ cc These two measurements are consistent with each other, though the central values 
are somewhat different. We shall take the averag ed value of 0.44 ± 0.09 for Dsj/D*tt°. The ratio 



DsTT+TT-/D*7r° is measured to be 0.14±0.04±0.02 by Belle |5l|. As for D*j/D*TT^it is found to be 
less than 0.22, 0.31 and 0.16, respectively, by BaBar 0, Belle Jsil and CLEO Q. Theoretically, 
the Ml transition Dgi turns out to be quite small Assuming that the I?si(2460) 

decay is saturated by D*it^, Dg^-, D*^ and DsTttt, we are led to 

0.53 < e(Z?,i(2460) ^ L>*7r°) < 0.68 . (3.44) 

This in turn implies B[B DL>si(2460)] = (1.6 ~ 4.6) x 10~^. As a result, the decay constant of 
£•^1(2460) is found to be 

/d.i (2460) -110 -190 MeV. (3.45) 

Our predictions fD*^ = 71 MeV and Jd^^ = 117 MeV with the latter being obtained from the 
relation 

fosi = cos Bs + /^3/2 sin 9s (3.46) 



with 6s ~ 7° inferred from the potential model are in agreement with experiment. 

Ideally, the neutral B decays B^ D**^n^ that receive only color-allowed contributions can 
be used to extract B D** transition form factors. Unfortunately, such decays have not yet been 



TABLE XIV: The predicted branching ratios for B^ D**^-k^ decays in the covariant light-front 
(CLE) and ISGW2 models. Since the decay constants of p-wave charmed mesons are not provided 
in the latter model, we employ the CLE decay constants and the ISGW2 form factors for the 
SGW2 results quoted below. Experimental results are taken from Behe BaBar and PDG 



18l [. The axial- vector meson mixing angle is taken to be = 12° 38] and the parameters ai^2 are 
given by ai = 1.07 and 02 = 0.27. 



Decay 


This work 


ISGW2 


Expt 


B- Z)(5(2308)°7r^ 


7.3 X 10^"^ 


4.8 X 10-^ 


(9.2 ±2.9) X 10-"^ \17j 


B- Di{2A27)°iT- 


4.6 X 10"^ 


9.4 X 10-^ 


(7.5 ± 1.7) X 10"^ [17] 


B- D[{2i20)°Tr- 


1.1 X 10-3 


8.2 X 10""^ 


(9.3 ± 1.4) X 10-^ [17, 54] 








(1.5 ±0.6) X 10-3 [18] 


B- D^{2m)^7r- 


1.0 X 10^3 


5.7 X 10"^ 


(7.4 ±0.8) X 10-^ [17, 54] 
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measured. Nevertheless, the decays B 



D**^iT that receive both contributions from color- 



allowed and color-suppressed diagrams provide a nice ground for testing the B D** form factors 
and the D** decay constants. Following jS^], we show the predicted branch ing ratios in Table [XTVI 



The experimental results are taken from Belle 17 1, BaBar [s^ and PDG 
we combine the Belle and BaBar measurements 



151. For B- 



B{B- 
B{B- 



^Df^~)B{Df 
Df^-)B{Df - 



^ D+TT~) 

D*+Ti-) 



(3.1 ±0.4) X 10-^ 
(1.8 ±0.4) X 10~^ 



to arrive at 



B{B- DfTT~)B{D^" D+TT-,D*+^-) = (4.9 ± 0.6) x 10" 



1*0 



(3.47) 



(3.48) 



Using B{D2^ — > D^TT ,D*^tt ) = 2/3 following from the assumption that the width is sat- 
urated by and D*tt, we are led to B{B- Dfix-) = (7.4 ± 0.8) x 10~^. We see from the 
Table [XTVI that the agreement between theory and experiment is generally good. In particular, the 
suppression of the D^'k~ production relative to the D^'7r~ one (the branching ratio for the latter 
being (5.3 ± 0.5) x 10^'^ UM) clearly indicates a smaller B Dq form factor relative to B ^ D 
one. For comparison, we also show the ISGW2 predictions in the same Table. Since the decay 
constants of p-wave charmed mesons are not provided in the ISGW2 model, we employ the decay 
constants in this work and the ISGW2 form factors to obtain the ISGW2 results quoted in Table 
IXIVI Predictions in the other models are summarized in (s^. 

Since the tensor meson cannot be produced from the V — A current, the decay B^ D2^tt~ can 



be used to determine the form factor combination r]{q^) = k(q'^) ± b+{q'^){m 



m 



' + b-{q^)q^ at 



The measured rate implies that r]{rn'^) = 0.43±0.02 , to be compared with the predictions 



of 0.52 and 0.38 in the covariant LF and ISGW2 models, respectively. 
It is worth mentioning that the ratio 

B{B- Dl{2mfTT') 



R 



B{B" D[{2420)0tt' 



(3.49) 



is measured to be 0.80±0.07±0.16 by BaBar |5J], 0.77±0.15 by Belle |l7| andl.8±0.8by CLEO |55|. 
The early prediction by Neubert jsfl yields a value of 0.35. The predictions of i? = 0.91 in the 
covariant light-front model and 0.67 in the ISGW2 model are in accordance with the data. 



IV. HEAVY QUARK LIMIT 

In the heavy quark limit, heavy quark symmetry (HQS) provides model-independent con- 
straints on the decay constants and form factors. For example, pseudoscalar and vector mesons 
would have the same decay constants and all the heavy-to-heavy mesonic decay form factors are 
reduced to some universal Isgur-Wise functions. Therefore, it is important to study the heavy 
quark limit behavior of these physical quantities to check the consistency of calculations. Since 
the analysis of heavy hadron structures and their dynamics in the infinite quark mass limit has 
been tremendously simplified by heavy quark symmetry and heavy quark effective theory (HQET) 
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developed from QCD in terms of l/mg expansion 53], it would be much simpler to study the 
decay constants and form factors directly within the framework of a covariant light-front model of 
heavy mesons fully based on HQS and HQET. Indeed, we have constructed such a model in 0] 
which can be viewed as the heavy quark limit of the covariant light-front approach discussed in Sec. 
II. We shall show explicitly that the decay constants and form factors obtained in the covariant 
light-front model and then extended to the heavy quark limit do agree with those derived directly 
in the light-front model based on HQET. 

Before proceeding, it is_ worth making a few remarks: (i) Just as in the conventional light-front 
model, it is assumed in 9] that the valance quarks of the meson are on their mass shell in the 
covariant light-front model based on HQET. However, this is not in contradiction to the covariant 
light-front approach discussed in Sec. II. As stressed before, the antiquark is on its mass shell after 
p~ integration in the covariant light-front calculation. Moreover, the off shellness of the heavy 
quark vanishes in the strict heavy quark limit. Therefore, the calculation based on the light-front 
model in 0] is covariant. (ii) Since the heavy quark-pair creation is forbidden in the mq oo 
limit, the Z-graph is no longer a problem in the reference frame where g"*" > 0. This allows us to 
compute the Isgur-Wise functions directly in the timelike region. 

In this work, we will adopt two different approaches to elaborate on the heavy quark limit 
behavior of physical quantities: one from top to bottom and the other from bottom to top. In 
the top-to-bottom approach, we will derive the decay constants and form factors in the covariant 
light-front model within HQET and obtain model-independent HQS relations. In the bottom-to- 
top approach, we study the heavy quark limit behavior of the decay constants and transition form 
factors of heavy mesons obtained in Sees. II and HI and show that they do match the covariant 
model results based on HQET. 



A. Heavy quark symmetry relations 

In the infinite quark mass limit, the decay constants of heavy mesons must satisfy the HQS 
relations given by Eq. 1)2. 7() . while all the heavy-to- heavy mesonic decay form factors are reduced 
to three universal Isgur-Wise (IW) functions, for s-wave to s-wave and ti/2 as well as for 
s-wave top-wave transitions. Specifically, B D,D* form factors are related to the IW function 
e(^) by El 

1 



1 



{niB +mD)f+{q ) + (m^ - mD)f-{q 



■- -2^mBmD* g{q^) 
a-{q')\ 



and obey two additional HQS relations 

a+{q^) + a-{q^) = 0, (mg - ■mD)f+{q^) + {mB + mD)f-{q^) = 0, 



(4.1) 



(4.2) 
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where lo 



■ m 



quark limit are related to Ti/2('^) via j25j 
n/2(w) 



q'^)/{2mBmn(*))- The B Dq and B D^'^ form factors in the heavy 



1 



{niB - mD*)u+{q^) + (ms + mDi)u-{q^) 



2^ /mBm^i/2 uj - 1 



^mBm^i/2 qi/2{q 



mBm^i/2 



while B 



U\l'^ and B 



T3/2(^^) 



-(cf(.^)-ci/^(.^)), 
transition form factors are related to ^"3/2(1^) by 
2 4/2(g^) 



(4.3) 



mBm^3/2 UJ 



\ 



2m| cf 



1 

3/2/9 



2m\ 



/ mB k{q ) 
3m £)* 1 + w 



uj-2 
2V2 



l+UJ 



^mBm^3/2 q3/2{q^) 



(4.4) 



and subject to the HQS relations 

b^{q') + b^{q')=0, c'f{q')+c]/\q')=0, 

{mB + mD*)u+{q'^) + {tub - ■mD*)u-{q^) = 0. (4.5) 

In terms of the dimensionless form factors defined in (|3.2|1 and (|3.6|) . Eqs. (|4.1|) and (|4.2j) can 

be recast to 



2^mBmD pBD^^2^ 
mB + rriD ^ 



2^mBmD F^^iq'^ 



niB + rriD 



1 



2^mBmD* yBD*^^2^ 
rriB + mo* 



mB+rriD* ° 



^2 [q ) 



mB + rriD* " ' uiB + mD 
Likewise, Eqs. 1)4. 3p and 1)4. 5() can be rewritten as 



(ms+mj3«)^ 



(4.6) 



Tl/2{^) 



^mBmB'^^BDl^^2^ 



iriB - niD* 



mB - mo* 



msm 1/2 



mB'm^i/2 



TUB - m 1/2 

^1 



niB - m^i/2 



TUB - m 1/2 

^1 



mBm j^i/2 



uiB - m^i/2 



,BD 



1/2 



1 



(4.7) 
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In next subsections we will derive the above HQS relations for form factors and IW functions using 
the covariant light-front model based on HQET. 

wwyrm 1/2 

We see from Table IVTni that HQS relations (|4.5|) for form factors b± and cj. are respected even 
for finite heavy quark masses. From the numerical results of Ti/2('^ = 1) = 0-31 and T^/2i^ = 1) = 
0.61 to be presented below in Sec. VI. F, one can check the HQS relations (|4.7() and 1)4. 4() at the 
zero-recoil point. It turns out that, among the fourteen B — > D** form factors, the covariant light- 

p;„i/2(3/2) BD^^^ BD^^^ 

front model predictions for A i ,Vq ^ , F2 ^ ,h,b+,b- are in good agreement with those 

3/2 3/2 

in the heavy quark limit, while the agreement is fair for and c_ . However, the predictions for 

BD* BZ)V2{3/2) 

-^10°'^! ^ ^ ^^^^ recoil show a large deviation from the HQS expectation. Indeed, 

Eqs. (|4.7() and ()4.4() indicate that except for ° , these form factors should approach to zero when 
q'^ reaches its maximum value, a feature not borne out in the covariant light-front calculations for 
finite quark masses. This may signal that Aqcd /w-Q corrections are particularly important in this 
case. Phenomenologically, it is thus dangerous to determine all the form factors directly from the 
IW functions and HQS relations since l/mg corrections may play an essential role for some of 
them and the choice of the /? parameters for s-wave and p-wave wave functions will affect the IW 
functions. 



B. Covariant light-front model within HQET 

To begin with, we rescale the bound state of a heavy meson by {P^ ^ Ph±i J, Jz) = 
\/Mh \ H{v, J,Jz))- It is well known that in the heavy quark limit, the heavy quark propagator can 
be replaced by 

i ^ i(l+ ^) ^ 

fig — tuq + ie 2v ■ k + ie^ 

where pq = rnqv + k and k is the residual momentum of the heavy quark. One can then redo all the 
calculations in Sec. II by using the above propagator for q[ and and perform the contour integral 
as before. Since the contour integral forces the antiquark to be on its mass shell, it is equivalent 
to using the so-called on-shell Feynman rules 3 in calculations. The zero mode effect arises from 
the = p'{^ = region and it can be interpreted as virtual pair creation processes . In the 
infinite quark mass limit, both quarks are close to their mass shell and far from the p'i~ = pl'^ = 
region. Consequently, the pair creation is forbidden and the zero mode contribution vanishes in the 
heavy quark limit. Hence, we do not need to stick to the = frame and are able to study form 
factors directly in the timelike region. 

To extract the on-shell Feynman rules, we use the calculation of the pseudoscalar meson anni- 
hilation [c.f. Fig. 1(a)] as an illustration. By virtue of Eq. ()4.8|) . the matrix element of Eq. (|2.8() 
can be rewritten as 

AT f H' 

where we have used p'^ = pq, p2 = Pq, m2 = ruq, N2 = Pg — + ie and k = —pq + {Mh — 'mq)v 
from 4-momentum conservation. As in Sec. II, we need to perform the contour integral by closing 
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the upper complex -p^ -plane, or equivalently, the lower complex -plane. The integration forces 
= vr?q and consequently, 



(4.10) 



Since X2 is of order Aqcd/?71q in the heavy quark limit, it is useful to define X = mQX2, which is 
of order Aqcd even if mq — > cxd. For on-shell Pq we have p~ = + m'^)/p'^ and 

1 

2X' 



v-p, = :^[pl+ml+X^). (4.11) 



It is then straightforward to obtain 



- , ^ — ^'(X.pj), (4.12) 

with the aid of Eqs. 1)2. 2p . 1)2.11(1 and the replacements 

Mo ^2{v -pq + mq)mQ , 

^i^,pl) - \l^m,Pl)- (4.13) 

An important feature of the covariant model is the requirement that the light-front wave function 
must be a function of v ■ Pq 0|: 

HX,pl)^^vpq). (4.14) 

As we will see later, the widely used Gaussian-type wave functions have such a structure in the 
heavy quark limit, while the BSW wave function does not have one. The normalization condition 
of $(u • Pq) can be recast in a covariant form: 



(2vr) 



^{27T)6{p'q-m'q)mvpq)\' = l, (4.15) 



or 



/roo 
X J 2(2^) 
Putting everything together we have 

{2n) 



f ''P^-mX,pl)\' = l- (4.16) 



(0|^^|P(^)) = -i^^ld%{-2^t)5{pl-ml) 



X -^^=^^$(^.p,)Tr[7^75(l+^)75(-^, + m,)]. (4.17) 
In practice, we can use the following on-shell Feynman rules to obtain the above and other ampli- 



tudes. The diagrammatic rule is given as follows 

(i) The heavy meson bound state in the heavy quark limit gives a vertex (wave function) as 
follows: 



-LJ ^{vpq)iTH, (4.18) 

V A^c V V -Pq + mq 
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'^*ivpg)iTH, (4.19) 



with Th = 7'^rJ:^7'^. For p-wave mesons, we denote the covariant wave function by ^p{v ■ Pg) 
(ii) The internal hne attached to the bound state gives an on-mass-shell propagator, 

'1+ 

(for heavy quarks), (4.20) 



« : i{— + mq) (for light antiquarks), (4.21) 

where v"^ = 1 and Pq = m^. 

(iii) For the internal antiquark line attached to the bound state, sum over helicity and integrate 
the internal momentum using 

d'^Pq 



Nr. 



-X-27ri)6{pi - mi) 



(4.22) 



where the delta function comes from the on-mass-shell condition and Nc comes from the color 
summation. 

(iv) For all other lines and vertices that do not attach to the bound states, the diagrammatic 
rules are the same as the Feynman rules in the conventional field theory. 

These are the basic rules for the subsequent evaluations in the covariant model. The vertex Th 
for the incoming heavy meson can be read from Table I or from Eqs. (|A11|) and HA19|) by applying 
Eq. (|4.13|) . Hence, the vertex functions in the heavy quark limit have the expressions: 



'Pi 



'P 



3r 



1 

V2 



-«75 
-^{V -Pq + mq 
{V-Pq^mq){4 



(4.23) 



75 



In terms of the and P^'' states, the relevant vertex functions read 



,3/2 



Pi 



1/2 



PI 



3/2 _ 



■^{V -Pq^mq) ^75 



(4.24) 



C. Decay constants 

In the infinite quark mass limit, the decay constants are defined by 

mi^l^K\P{v)) = iFpv^' , {0\qj^K\P*{v,e)) = Fye'' , 
mrK\S{v,e)) = Fsv^ mr-i^K\A^'\v,e)) = F^v2e^ 
mi^l^hM'^^v.E)) = F^3/2e^ (4.25) 
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where the decay constant Fh is related to the usual one fn by Fh = ^/MiLIh- Note that the 
tensor meson cannot be created from the V — A current. HQS demands that |2ll. 



Fv = Fp, = Fs, F^3/2 = 0. (4.26) 

Using the Feynman rules shown above, it is ready to evaluate the one-body matrix elements for 
heavy scalar and axial- vector mesons: 

{0\qj''^5K\A'/\v,e)) = -l=Tl■{^^^J5^ ^l^Mi} , (4.27) 

where 



Ml = VN, I ^(27r)<^(p^ -m^)-^=^|=|^(m,- ^^)iv . + m,) . (4.28) 



Letting Aii = ai + bi ^, we obtain 



ai = Vn'c I 7— ^(27r)5(pg - ml)<^>p{v ■ Pg) rUg ^/v-pq + ruq, 



{2n) 



d Pq 



Thus, 



bi = -VNcJ j^i2Tr)5ip'g - m^)$p(i; • pg)iv ■ pg)^v ■ Pg + nig . (4.29) 
2 



Nc f d'^Pg 




3 J (27r) 



^{27r)6{pg - mg)<^>p{v ■ Pg)^v ■ Pg + nig {v ■ Pg - nig) 



3 J 2{2TTfX '^-^'V 2X 2X 



Likewise, for the axial-vector meson 



{0\qri5hM'^\v,e)) = -l^Tt{rj,^[?!{^^ + v^) - Se^h^M^}, (4.31) 



V6 I' 2 

with 



M^ = ^el 7;^(27r)5(p2 _ m^) J\ ^'^^ {nig- ^,)p°. (4.32) 



P^A2.)5{pl-nil)4M=^. 
(2vr)4^ > i>^vpg + nig 

The general expression of is 

M'^ = a^v'^ + 627" + C2 ^v" + d2 ^7" • (4.33) 

Since e ■ v = Q and the contraction of 7^ with the spin 3/2 field vanishes, namely, 

(1+ ^) [^(7a + v^)- 3e„]7" = 0, (4.34) 

we are led to -^^3/2 = in the heavy quark limit. 

For completeness, the decay constants of the s-wave heavy mesons are included here 

Fp = Fv 
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We now show that the decay constants obtained in the covariant hght-front model in Sec. II do 
respect the heavy quark symmetry relations given in ()2.7I) or (|4.26|) in the infinite quark mass hmit 
and have expressions in agreement with Eqs. 1)4. 3U() and 1)4. 35() . To illustrate this, we consider the 
decay constants of pseudoscalar and vector mesons. In the m'l = mq oo hmit, Eqs. (|2.17j) and 
1)2. 22() are reduced to 

, , , /iVc f dXdPp^ 2^ X + rriq 




where niq = m2, X = mqx and use of Eq. (|4.13|) has been made. Since the wave function is even 
in a quantity defined in Eq. ()2.2() . it fohows that 



j dx(fp± 



■■Pz 



[rnqx + mq{l - x)]^ + p'j_ 

Therefore, fv = fp in the heavy quark limit. Moreover, ^/fnqfp is identical to Fp in Eq. I|4.35j) 
after applying the identity 



m.A) . + = / HX.pl) ^^''^'"f-^A^ (4.38) 



fohowing from Eq. (|4.37j) . 

Likewise, Eq. ()2.23|) for the decay constants of axial-vector mesons is reduced in the heavy quark 
limit to 

^''^ ^ 'VTi 2(2.) VX ^ix^,u,y^f/ '' 



2(2»)3v/X ^(.Y + m,)2+pi V X^- 

where use of Eq. 1)4. 37(1 has been applied for deriving the expression of /i^. By virtue of the 
identities 

dXd'^p^ %{X,pI) 



J{X + mqY+pl 



{pi- pi/2) 

+pi 

dXd-p^ %{X^p\) (440) 



and 



IX . /(X + m,)2 + p\ 

3 fdXd'ps. MX,pI) 2 
' P± 



V2J VX J^x + mqf+pl 
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following from the first equation of (|4.4UI) . one can show that /s^ = — -v/2/i^ and hence /y^3/2 = 
and ^/rriQf ^1/2 = = Fs in the mq oo limit. 

D. Isgur-Wise functions 

It is well known that the s-wave to s-wave meson transition in the heavy quark limit is governed 
by a single universal IW function ^(to) j2]|. Likewise, there exist two universal functions Ti/2{(^) 



and T3/2(w) describing ground-state s-wave to p-wave transitions |25|. Since the IW function ^ has 
been discussed in detail in 0|, we will focus on the other two IW functions T1/2 

1/2 

Let us first consider the function T1/2, which can be extracted from the B — > Dq or B ^ D^' 
transition 

{Dl{v')\K,Thl\B{v)) = -i-LTr|(^l^)r(^i±^)75-M3}, 

{D\'\v\e)\K,Thl\B{v)) = i-LTr|^75(^)r(i^)75-M3}, (4.42) 
where 7W3 is the transition matrix element for the light antiquark: 

-M3 = / [d%]{niq- ^q){v' ■ p, + mq) , (4.43) 
and we have introduced the short-hand notation 

. - ,4.44) 

^ ' \J {V ■ Pq + mq){v Pq + mq) 



The structure of A^3 dictated by Lorentz invariance has the form JS^ 

7W3 = 03 + 63 ^ + C3 / + ci3 ^ / . (4.45) 

This covariant decomposition allows us to easily determine the coefficients 03,63,03,^3 with the 
results: 



03 = y" [d'^Pq] mq{v' -Pq + mq), 



[<iVj W ■ P, + mM (^i±^ + <^i^^ \ . (4.46) 



C3 = - y [d%] {v' ■ Pq + mq)- 

ds = 0, 



ij {V + V') ■ Pq {V - V') ■ Pq 



1 + UJ 1 — UJ 



where uj = v ■ v' . 



1/2 

Then B Dq and B ^ transitions are simplified to 



{D*iv')\K,rhl\Biv)) = -^2ri/2(^)TY|(i^)r(i^)75~ 
{Dl/\v',e)\K,rh',\B{v)) = i2n/2(^)Tr|/75(^)r(ii^)75}, (4.47) 
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with 



Ti/2(w) = ^^(03-63 + 03-^3) 



Since 



= ^ / [d%] {v' ■ P, + m,) (m, + ■ (4.48) 



v-p, = ^{pl + ml + X% v'-pg = ^{pI + ml + X% (4.49) 



the IW function T1/2 can be exphcitly expressed as 



= /life 



V pi + {mq + Xy 
where z = X' /X and it can be related to v ■ v' by 



2V3J 2(27r)3X2 v^(l-z) 

pi + (mg + X)(mq-zX) , (4.50) 



IpI + {nig + zX)2 



± Vu;2 - 1 , z+ = —, (4.51) 

with the + (— ) sign corresponding to greater (less) than v'^. Note that greater (less) than v'^ 
corresponds the daughter meson recoiling in the negative (positive) z direction in the rest frame 
of the parent meson. In other words, after setting v± = v'j_ = 0, the daughter meson recoiling in 
the positive and negative z directions arc the only two possible choices of Lorentz frames. It is 
easily seen that Ti/2{^) remains the same under the replacement of z ^ 1/z. This indicates that 
the Isgur-Wise function thus obtained is independent of the recoiling direction, namely, it is truly 
Lorentz invariant. 

Next consider the 5 — Z)| or S — D"^^^ transition to extract the second universal function 

{D*2iv',s)\K.Thl\B{v)) = -iTr 1^8^^^ (^) T TsA^f | , (4.52) 

{D',/\v',e)\K,Thl\Biv)) = --!=Tv{[(-7a + 0^+3£a]75(^)r(^)75A^?|, 
where 

M'i = l[d%]{mg- fig)p^. (4.53) 

Its most general expression is 

M^a = a^Va + hv'a + C4 '^Va + /va + 64 J^v'^ + /v'^ + g^Ja + ^4 J^la + ^4 /Tq- (4.54) 

Although only terms proportional to 04, C4 and ^4 will contribute to B ^ D2 and B d\^'^ 
transitions after contracting with the vertex of the spin 3/2 particles, all the terms in M.ia have to 
be retained in order to project out the coefficients. With M.Aa contracting with v°^, v"^ and 7° we 
find the following equations: 



04 + 64W = y \d^Pq\ rrigV ■ pq, 
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d^Pg] ml, 



d'pg] 



'd'^Pg] 



a^uj + &4 = 

C4 + 2diOJ + /4 + 4^4 = 

C4 + d^u: + (74 = 

d4 + /4a; = 

C4W + 6^4 = 
^4^*^ + /4 + 54 = 

and 64 = (i4, h4 = h'^ = 0. Solving the above equations yields 

«4 = ^ /[A,](A+ + A_), 



V ■ Pq ( {V + v') ■ Pg {V - v') ■ Pg 



l+OJ 



+ 



l-OJ 



4 , f • / (f + w') • Pq {V - V') ■ Pg 



l+LO 



l-UJ 



V' -pg ({V + V')-Pq {V-V')-Pg 



l+UJ 



+ 



l-UJ 



4 -Pg f{V+v')-Pg {V-V')-Pg 



d'pg] 



l-UJ 



hi 

C4 

di 

h 



nin 



[d%](A+-A_), 



-\j[d%] ((A+ + A_)2- 
\j[d%] ((Ai-A^_) + 
-\j[d%] (^(A+-A_)2- 



{l+u:){l-Lo) 

{i + uj){i-ujy 
94 



and 



with 



(l+w)(l-w) 
9^ = -\J id%] (rnl - 1(1 + ^)A^ - ^(1 - lo)XI 



A_ 



{v + i;') • Pg 



A_ 



(U - V') - Pg 



(4.55) 



(4.56) 



(4.57) 



(4.58) 



I + UJ l-UJ 

Since only v^, and ^'v^ terms in M4Q, survive after contracting with the vertex of D2 and 
Z?^^^ particles, the matrix elements (|4.52() are simplified to 



{Dl{v',e)\K,Th':,\B{v)) = -^^/3T3/2(^)^./3^;^Tr|7"(^)^(^i± 
{D',/\v',e)\K,rhi\B{v)) = -^ry.iu) 

X Tt S^y{l + uj) + 3e ■ v]ji 



75 



(4.59) 



fl+ f\^fl+ 



V 2 



V 2 



75 



with 



^3/2 ( 



1 

7i 
1 

2\/3 



(04 — C4 — ^4) 



(4.60) 



(A+ + A_)(mg + A+) + 



A^(l + a;) + A^(l-a;)-m^ 
2(1 +u;) 
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A more explicit expression of T3/2 reads 
, , If dX(Pp± 

1 



[2(l-w)(l+a;)2 
+ 3{v pgf - {I - uj^)m: 



2VI^;{zX,pI)^X,pI) 
[p\ + (mg + X)2][p2_ + + zXy^ 

{1-2UJ){V' ■Pg){2V + v')-Pg 



+ 



1 



[v -pq- u;{v' ■Pq)]mg 



1_^2 

After some manipulation, we obtain a simple relation between r3/2 and ti/2- 



UJ 



/[p'_2 + ^rng + X)2][p'_2 + (mg + zX)^] 
Finally we include the usual IW function S,{uj) for the sake of completeness 

dX(fpx. 2y/l 



, ^*{zX,p\)^{X,p\) 
p\ + {mg + X){mg + zX) 



'\pI + {mg + XY]\p\ + (mg + zXY\ 
and the relevant matrix elements are given by 

{D{v')\hl,Thi\B{v)) = ^(a;)Tr{75(^)r(^)75} , 



(4.61) 



(4.62) 



(4.63) 



(4.64) 



E. Form factors in the heavy quark limit 

From Eqs. (|4.64l) . (|4.47|) and (|4.59|) we obtain the matrix elements of B ^ D, D* ,D** transitions 



in the heavy quark limit 
{D{v')\V, 
{D*{v\e)\V, 
{D*{v\e)\A^ 
{Dl(v')\A, 
{d\'\v',e)\V, 
{D\/\v\e)\A, 
{Dl'\v',e)\V, 

{Dl'\v\e)\A, 

{Dl{v\e)\V, 
{Dl{v\e)\A, 



B{v 
B{v 
B{v 
B{v 
B{v 
B{v 
B{v 

B{v 

B{v 
B{v 



i2Tii2{u)){v - v')/,, 
-i2Ti/2{uj)Ul-u;)e*^ + {e* •v)v'^ 



-2ri/2(a;)6^,„^e*^-'"-^ 



V V'' 



(4.65) 



) = i^Ty2ii^){{l-co')s;-ie*-v)[3v, + i2-co)v'^]}., 

1 



7-3/2('^)(l + i0)ef,uaf3£*''v'"v'^, 

V3T-^/2{uj)e^,yape*'''^Vjv"'v'^ , 
-iv^r3/2(u;){(l + u;)e;X " 4;3^"^^^^m}- 
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It is easily seen that the B D** matrix elements of weak currents vanish at the zero recoil point 
uj = 1 owing to the orthogonality of the wave functions of B and D** . Setting pB = rnsv and 
Pd = rriDv' ■ •, etc. in Eqs. and 1)3. 5|) and comparing with (|4.65|) yields all the form-factor 

HQS relations given in Sec. IV.A. 

We are ready to check the heavy quark limit behavior of form factors to see if they satisfy the 
HQS constraints. Consider the form factor °{q'^) = —u^^q^) first. Let X2 = x,xi = l — x,m2 = 



niqjX = miyX,X' = nicX, it follows from Eq. (|3.2U|) that 

n+{g') = ^ldxdy^^^{2iX + m,)X'-2m,X-2mlil-xf-xV 

+2mgx{X - X') - 2pl + 2xp_^ • g_L }, (4.66) 

where use of Eq. ()2.2|) and ~ P'± ~ ^Q-L been used. In the heavy quark limit 
X ~ ©(Aqcd/wiq) 0, we have 

Substituting the replacements 



VMI) - sj^HX.p'l), V,(x,v'l) - ^|^^p{X',p'l) (4.68) 
valid in the infinite quark mass limit and noting that z = X' / X = mc/nrn,, we arrive at 



and hence 



2^7 2(27r)3X2 z P 
/|±K±|f[,«,(„, + ^„„.,_,^,]. ,4.60) 



(4.70) 

ruB - mo* 



Likewise, it is easily shown that 

^'™^-Fr(g^)-C(^). (4.71) 



rriB + mo 

In order to demonstrate that the B D* form factors are related to the IW function S,{uj), we 
need to apply the identity (|3.36jl which has the expression 

dxd^^^^i^ ( {q . P)^^ - -ip'l + ml- X^)] = (4.72) 

in the mq oo limit. This identity allows us to integrate out the (p^ • qi_)/q^ term. Then the 
form factor g{q^) that reads [see Eq. (|B4I) ] 
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is reduced under the heavy quark hmit to 



where use of Eqs. H2.11() . (|3.36p and ()4.68p has been made. Comparing with Eq. H4.63() it is evident 
that the heavy quark hmit of g{q'^) has the same expression as ^(w) apart from a mass factor. 
Therefore, we arrive at'' 

^V^^V^^\q') = -2.M^9iq') - Ciu^)- (4.77) 
niB + rriD* 

We next turn to the form factors qi/2 and see if they are related to the IW functions 

Ti/2 and T3/2, respectively. We first study the heavy quark limit behavior of q^^ and q^^. It follows 
from Eq. iHO)!) that 



msm^^/.q (q ) - ^'^Bm^^f2J^ j dxd P±^J^ [P± + 



1 f dXd'p'^ ^,, 

2 J 2(2^)S 



and 



[p'l + {m, + Xnp'l + {mg + zX)^] 

1 f dxdy^ %{zx,p'i)^{x,p'i) 



(4.78) 



p'l + (mg + zXy r ,2 , , 



Since 



. (4.79) 

/[p'2 + (^^ + X)2][p'_2 + (m, + zXy] 

^ If tlie Taylor expansion of hy/N[' is performed to take care of the p'j_ ■ q± term in the integrand of g{q^), 
it turns out that the heavy quark limit of the B D* form factors will be related to the IW function 

7 2(27r)^ ^[p'J^ + {niq + X)^][p'l + (m, + X'^] 

where 

/V" / H h" \ 

^v-Trr{-rJi2^) • (4-76) 

This function C('^) first obtained in was found numerically identical to ^(w), as it should be. However, 
one has to appeal to the identity H4.72|l in order to prove this equivalence analytically. 
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following from Eq. H2.6|) . we obtain 

^msm^i/a qi/2{q^) Ti/2{uj) (4.81) 

and 

- YTu}\I^^^J^T ^3/2(9^) ^ '^3/2(t^), (4.82) 

as promised before. 

All other HQS relations in Eqs. ()4.1() . 1)4. 3p and 1)4. 4j) can be proved in the same manner except 
for the B D2 form factors h, k, b^, b- for which we are not able to show at present that they are 
related to T3/2{^) in the heavy quark limit. Perhaps one needs some identities in and those 
derived in Appendix B to verify the HQS relations between B D2 form factors and T3/2{^)- This 
remains to be investigated. 



Numerical results for IW functions and discussion 



Covariance requires that light-front wave functions be a function v -pg. Currently, there exist 
several phenomenological light-front wave functions commonly utilized in the literature. There are 
several popular phenomenological light-front wave functions that have been employed to describe 
various hadronic structures in the literature. Two of them, the Bauer-Stech-Wirbel (BSW) wave 
function ^BSwi^jP^) [i^] and the Gaussian-type wave function ^ci^^P^) [^B; have been widely 
used in the study of heavy mesons. In the heavy quark limit, we denote these wave functions as 
follows: 




X exp 



(4.83) 



where $g(-'^iPj_) is the heavy quark limit expression of the Gaussian- type wave function given in 
Eq. (|2.12|) . For p-wave heavy mesons, the wave functions are 



BSW 



{X,pi), <^^{X,pi) 



ci>G(x,p1). 



(4.84) 



As pointed out in 9], not all the phenomenological light-front wave functions have the covariant 
property. We found that the Gaussian wave function and the invariant-mass wave function can be 
reexpressed as a pure function of v ■ pg. The wave function can be rewritten in terms of v ■ pg. 





m„ 



■ Pg exp 



V-PqY 



m„ 



(4.85) 
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FIG. 6: The Isgur-Wise functions ^, rx/2 and t-^i2 as a function of 



ijj. 



Therefore this wave function preserves the Lorentz covariance of Eqs. (|4.54)) . (|4.45|) and (|4.28j) . 
This also can be examined by a numerical check of the covariant condition 



cD(X,pi)X 



(4.86) 



2(27r)3X ' 7 2(27r)3X X 

which is satisfied if <I>(X, p^) is a function oiv -"Pq. However, very surprisingly, the commonly used 
BSW wave function cannot be recast as a pure function of u • pg. Hence the BSW wave function 
breaks the Lorentz covariance. Indeed, we have already found previously that there is some 
inconsistent problem by using the BSW wave function to calculate various transition form factors. 
Now we can understand why the BSW wave function gives such results inconsistent with HQS 
found in 0, E^- Hence, by demanding relativistic covariance, we can rule out certain types of 
heavy meson light-front wave functions. 

To perform numerical calculations of the decay constants and IW functions in the heavy quark 
limit, we follow 0| to use the input = 250 MeV and = 180 MeV to fix the parameter /3oo to 
be 0.49. For decay constants we then obtain 

Fp = Fv = 413 MeV^/^ Fs = F41/2 = 399 MeV^/^ 



(4.87) 



3 /2 

The decay constant of the heavy meson vanishes in the infinite quark mass limit. We see that 
the decay constants of ground-state s-wave mesons and low-lying p-wave mesons are similar in the 
heavy quark limit. 

The IW functions (|I3n|) and (|^^^ can be fitted nicely to the form 

,2 



and it is found that (see Fig. 

e(^) ^ 

7"3/2('^) ■ 



/(l)[l-p2(^_l) + ^(^ 



11 



(4.88) 



1 - 1.22(w - 1) + 0.85(w - if, 
0.31 (1 - 1.18(a; - 1) + 0.87(w - 1) 

1 - 1.73(0; - 1) + 1.46(a; - 1) 



0.61 ( 



(4.89) 
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TABLE XV: The Isgur-Wise functions T3/2 at zero recoil and their slope parameters. The 



numerical results for 



63 



6J] denoted by "*" are quoted from |6| 



Tl/2(1) 




P1I2 


-^3/2(1) 


2 

^3/2 


Ref. 


0.31 


1.18 


0.61 


1.73 


This work 


0.06 


0.73 


0.52 


1.45 


[61j* 


0.09 


1.1 


0.28 


0.9 


m 


0.13 


0.57 


0.43 


1.39 


[63j* 


0.22 


0.83 


0.54 


1.50 


[64]* 


0.34 


1.08 


0.59 


1.76 


[19]* 


0.35 ± 0.08 


2.5 ± 1.0 






[65] 


0.41 ± 0.04 


1.30 ±0.23 


0.66 ± 0.02 


1.93 ±0.16 


m 






0.74 ±0.15 


0.90 ±0.05 


m 



where we have used the same fioo parameter for both wave functions $ and It should be 
stressed that unlike Ti/2(1) and T3/2(1), the normalization ^(1) = 1 at the zero recoil point is 
a model-independent consequence; that is, it is independent of the structure of wave functions. 
In Table IXVl we have compared this work for the IW functions Ti/2(^) and t^/2{'^) with other 
model calculations. It turns out that our results are similar to that obtained in the ISGW model 
[1^ (numerical results for the latter being quoted from Our result = 1.22 for the slope 

parameter is consistent with the current world average of 1.44 ± 0.14 extracted from exclusive 
semileptoic B decays 

It is interesting to notice that there is a sum rule derived by Uraltsev j^] 

(n)/-^l') 1 



E 



4' 



(4.90) 



where n stands for radial excitations. This sum rule clearly implies that |t3/2(1)| ^ I ''"1/2 (1)1- Our 
results indicate that this sum rule is slightly over-saturated even by n = p-wave states. Another 
sum rule due to Bjorken 



221 1 reads 
1 



+ EI 



where is the slope of the IW function ^(w 
to 



)P+2^|T<;^{1)|', (4.91) 

n 

Combined with the Uraltsev sum rule 1)4. 90(1 leads 



An) 

h/2 



(1) 



(4.92) 



Note that while the Bjorken sum rule receives perturbative corrections 7]|, the Uraltsev sum rule 
does not (for a recent study, see [T^]). 



48 



V. CONCLUSIONS 

In this work we have studied the decay constants and form factors of the ground-state s-wave 
and low-lying p-wave mesons within a covariant light-front approach. This formalism that preserves 
the Lorentz covariance in the light-front framework has been developed and applied successfully to 
describe various properties of pseudoscalar and vector mesons. One of our main goals is to extend 
this approach to the p-wave meson case. Our main results are as follows: 

• The main ingredients of the covariant light-front model, namely, the vertex functions, are 
explicitly worked out for both s-wave and p-wave mesons. 

• The decay constant of light scalar mesons is largely suppressed relative to that of the pseu- 
doscalar mesons and this suppression becomes less effective for heavy scalar resonances. The 
predicted decay constants I/d*,,! = 71 MeV and |/dsi{2460)I = MeV are consistent with 
the corresponding values of 47 ~ 73 MeV and 110 ~ 190 MeV inferred from the measurement 
of DD*Q and DDgi productions in B decays. 

• In the limit of SU(N)-flavor symmetry, the decay constants of the scalar meson and the ^Pi 
axial-vector meson are found to be vanished, as it should be. 

• The analytic expressions for P ^ S,A transition form factors can be obtained from that of 
P ^ P,V ones by some simple replacements. We have also worked out the form factors in 
P ^ T transitions. 

• The momentum dependence of the physical form factors is determined by first fitting the form 
factors obtained in the spacelike region to a 3-parameter function in and then analytically 
continuing them to the timelike region. Some of the V2(g^) form factors in P — ^ ^ transitions 
are fitted to a different 3-parameter form so that the fit parameters are stable within the 
chosen range. 

• Numerical results of the form factors for B{D) tt, p, ao(1450), ai(1260), 6i(1235), 02(1320), 
B{D) K, K*, K** and B ^ D, D*, D** transitions are presented in detail, where K** 
and D** denote generically p-wave strange and charmed mesons, respectively. 

• Comparison of this work with the ISGW2 model based on the nonrelativistic constituent 
quark picture is made for B D** transition form factors. In general, the form factors 

at small in both models agree within 40%. However, Fq and ^ {(^) have a 

very different behavior in these two models as (f increases. Relativistic effects are mild in 
B D** transitions but can manifest in heavy-to-light transitions at maximum recoil. For 
example, Vf^""^ (0) is found to be 0.13 in the covariant LF model, while it is as big as 1.01 in 
the ISGW2 model. 

• The decay amplitudes of B~ D**^7r~ involve the B D** form factors and D** decay 
constants. We have compared the model calculations with experiment and found a good 
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agreement. In particular, the suppression of the DqTt production relative to D tt one 
clearly indicates a smaller B ^ Dq form factor relative to the B ^ D one. 

• The heavy quark limit behavior of decay constants and form factors is examined and it is 
found that the requirement of heavy quark symmetry is satisfied. 

• Decay constants and form factors are also evaluated independently in a covariant light- 
front formulism within the framework of heavy quark effective theory. The resultant decay 
constants and form factors agree with those obtained from the covariant light-front model and 
then extended to the heavy quark limit. The universal Isgur-Wise functions ^{i^),Ti/2{'^) and 
Ts/2{i^) are obtained and a relation between T1/2 and is found. In the infinite quark mass 
limit, all the form factors are related to the Isgur-Wise functions. In addition to ^(1) = 1 
at zero recoil w = 1, it is found that ri/2(l) = 0.61, r3/2(l) = 0.31 and p"^ = 1.22 for the 
slope parameter of ^(w). The Bjorken and Uraltsev sum rules for the Isgur-Wise functions 
are fairly satisfied. 
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APPENDIX A: VERTEX FUNCTIONS IN THE CONVENTIONAL LIGHT- 
FRONT APPROACH 



In the conventional light-front approach, a meson bound state consisting of a quark giand an 
antiquark q2 with the total momentum P and spin J can be written as (see, for example j7|) 

|M(P,2^+l Lj, J,)) = J {d^pl}{d^P2} 2{27Tf6^P -pi- P2) 

X ^'ls{puP2,XiA2) ki(m,Ai)^2(P2,A2)>, (Al) 



where pi and p2 are the on-mass-shell light-front momenta, 



2 2 

P = [P ,P±) , P± = [P ,P ) , P = T — , (A2) 



and 

3 _ dp^(fpx. 

Ai)g-(P2, A2)) = 6i^(pi)4^(p2)|0), (A3) 
{&A'(P'),4(P)} = {^iA'(P ),4(P)} = 2(2^)=^ 6\p'-p) 5yx. 
In terms of the light-front relative momentum variables {x,p±) defined by 

pf = xiP'^, P2 = X2P^ ■, 2;i+X2 = l, 

Pl± = XiP^ +p_L, P2± = X2P± - P±, (A4) 

the momentum-space wave- function ^^5"" for a '^^'^^Lj meson can be expressed as 

(^i'P2,Ai,A2) = -^{LS;L,S,\LS;JJ,)RfX{x,p^) ^^llA^^Pi.), (A5) 

where ^LL^{x,p_[_) describes the momentum distribution of the constituent quarks in the bound 
state with the orbital angular momentum L, {LS; LzSz\LS; JJz) is the corresponding Clebsch- 
Gordan coefficient and Pf^x^ constructs a state of definite spin {S,Sz) out of light-front helicity 
(Ai,A2) eigenstates. Explicitly, 

-^AfA2(^'^'-L) = H (Ai|'^M(l-a^'l'-L>"^i)l'Si>(A2|'^M(a;,-p±,"^2)|s2) SiS2\^^; SSz\ , (A6) 

where \si) are the usual Pauli spinors, and TZm is the Melosh transformation operator 0, E^: 

UDiPi,s)u{pi,X) VD{Pi,s)v{pi,X) 



{s\'RM{x,p±,mi)\X) 



2mi 2mi 



m + XjMo + iasx -px xn 
{mi + XiMof + pj_ 



with tt(£)), a Dirac spinor in the light-front (instant) form, n = (0,0,1), a unit vector in the 
z-direction, and [cf. Eq. (|2.2I) ] 

M^,=^d±A + !!i±A. (A8) 

Xi X2 
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Note that uo{p,s) = u{p, X){X\TZ\j\s) and, consequently, the state \q{p, X)){X\TZ\.j\s) transforms 
Uke \q{p,s)) under rotation, i.e. its transformation does not depend on its momentum. 



SS^ 



In practice it is more convenient to use the covariant form for R^'^ 



\2 



RfXix,p^) 



V2 Mo (Mo 



uipi,Xi)iP + Mo)Tvip2,X2 



(A9) 



with 



Mo = ^M2 - (mi -m2)2, 

P =Pl +P2, 



e^(O) = 



P- 



:e^{±l) ■ P^, 0, e^{±l) 



el(±l) = T(l,±i)/V2 



1 

Mo 



P+ 



P^,Pi 



(AlO) 



For the pseudoscalar and vector mesons, we have 

Tp = 75 (pseudoscalar, 5 = 0), 

Tv = — ^{Sz) (vector, 5* = 1). 

It is instructive to derive the above expressions by using the relations 

ud{pi,si)ud{pi,si) 



(All) 



u{pi,Xi) 

V(P2,X2) 

Ud[PI, Si) 75Vd{P2, S2) 



u[Pi,X 



2mi 

VDiP2,S2)vDip2,S2) 

2m2 



V{P2,X2), 



UD{iei,p), Sl) ^'^^^ j5VD{{e2, -P},S2) 



^D(Pl,gl) '^^^^° (- ^iSz))vDiP2,S2) 



^ (ei + mi)(e2 + m2) ixl^cr2X*s2 
2(ei + mi)(e2 + ^2)1^ ^; siS2|^ \; 00^ , (A12) 

UD{{ei-,p)-,si) ^'^~^^ e{Sz) ■ 7VD{ie2, -p), S2) 



(ei + mi)(e2 + m2) ixl^e{Sz) ■ ^(72X1^ 



'2(ei + mi)(e2 + m2) 

where [cf. Eq. d^ ] 

Mo = 61+62, ei 



2(61 + mi)(62 + m2) 
Mo(Mo + mi + m2) 



i^i^;siS2\--;lSzn 



V2M0 



mf +p\+pj,, pz 



X2M0 ml+p\ 



(A13) 



2 2x2 Mo ' 

Xs is the usual Pauli spinor and we have used the usual properties, especially, the covariant one, of 
Dirac spinors. Applying equations of motion on spinors in Eq. ()A9|) leads to 



u{pi){P + Mq)^5v{p2) = (Mo + mi +m2)ti(pi)75u(p2), 

u{pi){P + Mq) ^v{p2) = u{pi)[{MQ + mi + m2) ^- e- {pi-p2)\v{p2) 



(A14) 
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and Rx^x2 reduced to a more familiar form ^3]. It is, however, more convenient to use the form 
shown in Eq. HA9|) when extending to the p-wave meson case. Two remarks are in order. First, 
Pi + P2 is not equal to the meson's four-momentum in the conventional LF approach as both the 
quark and antiquark are on-shell. On the contrary, the total four-momentum is conserved at each 
vertex in the covariant LF framework. Second, the longitudinal polarization 4- vector e'^(O) given 
above is not exactly the same as that of the vector meson and we have i ■ P = 0. We normalize the 
meson state as 



{M{P', f, Ji)\M{P, J, J,)) = 2{27rfP+6^{P' - P)5j,j5j,j^ , (A15) 

SO that 



/dx dF'pj_ 
2{2T:f 'p'l'L'S^^P±)vllAx,P±) = h',L (A16) 

Explicitly, we have 

(P00 = ^, 'PlL,=PL,Vp, (A17) 



where pl^=±i = T{p±x =t ip±y)/V^, PLz=o = Pz are proportional to the spherical harmonics 
in momentum space, and (/?, ipp are the distribution amplitudes of s-wave and p-wave mesons, 
respectively. For a Gaussian-like wave function as shown in Eq. (|2.12)) 7], one has ^pp = \/2/]3^(p. 
For p-wave mesons, it is straightforward to obtain 

{lS;L,S^\lS;JJ^)p,^Ri%{x,p^) = ~ ^ — - 

V2 Mo(Mo + mi+ 1712) 

xu{pu Xi)ip + Mo)r2s+ipjV{p2, As), (A18) 

with 

Tapo = 
Tip, = 

rap, = ^^(/r-^)^-e.i.J,„ 

Tap, = i^,r{-Kn, (A19) 

where K = {p2 —pi)/2, £fj,i,{m) = (11; m'm"|ll; 2m) e^(m')ei,(m"). Note that the polarization 
tensor of a tensor niGSon satisfiGS tliG relations! — 

and e^yP^ = = and that e^, e^y are 
identical to e^, respectively, for maximal transversepolarized states (m = itJ). The above 
expressions for ^Pi and ^Pq states are consistent with 0,01 and 69 1, respectively. 

The vertex functions shown in Table 1 and Eq. 1)2. 11() follow from the above explicit expressions 
for ^^ifg ■ For example, by taking e^y{—K'^) in place of in Eq. ()A14|) we obtain the '^Ps vertex 
in the form shown in Table 1. Note that there are an overall factor and sign to be determined. 
The overall factor, (M^ - M^).J^{^ [cf. Eq. dTTTl) ]. is fixed by comparing the pseudoscalar decay 
constant fp obtained in both covariant and conventional approaches [see Eqs. (|2.16|) and (|2.17|) ]. 
while the overall sign can be fixed by the HQS expectation for decay constants and form factors. 

1/2 Q 

For example, the sign of the P^ state relative to Pq is fixed by the HQS relation / 1/2 = fs- 
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An additional factor of i is assigned in Table 1 as in the usual Feynman rules to ensure that the 
corresponding operators are hermitian. For example, we have an i in front of 7^ but not 75, just 
like the usual QED and Yukawa vertices, respectively. Similarly, polarization vectors are decoupled 
from the vertex Feynman rules as usual. 



APPENDIX B: SOME USEFUL FORMULAS 



In this Appendix we first collect some formulas in relevant for the present work and then we 
proceed to summarize the formula for the product of four p'^s needed for the calculation in Sec. III. 



The explicit representation of the traces in Eqs. 1)3. 11(1 and H3.22() can be found in [l 
completeness we collect them in below: 



For 



~,pp 



+ 



//\2i 



2M'^ + N[- N'{ + 2N2 + 2{m[-m2)- {rn[ - m'l) 



+P^[q^ -N[-N'{-{m\-m'{) 



(Bl) 



and 



PV 



s7\ 



mi 



-2ie^,ap{pTP^imi 



+ p'^qf^im'i + m[ - 2m2) + q"P^m[} 



V 



J/2 



ml) 



mi) 



+2g^,u[m2{q^ -N[- N'{ - m'^ - mf ) - m'^iM"'^ - N'{ - N2 - m'[ 

— m"(M''^ — — — — m'l) — 2m'^m'/m2| 
+8p'ifj,Pi,yim2 - m'l) - 2{P^qy + q^Py + 2q^q,y)m'i + 2p'i^Py{m'i 
+Vi/i9i^(3m'i - m'{ - 2m2) + 2P^p'i^{m'i + m'{) + 2q^p'i^{?,m'i + m'{ - 2m2) 

+ ^(V. - 3(?. - P,){2p'i^[M'^ + M"2 - g2 _ 2N2 + 2{m'i - m2){m'i + ma)] 

- 2M'2 + N'i- N'; + 2N2 - {mi + m'/)^ + 2{m'i - ma)^] 
+P^[q^ -N'l- N'{ - {m'l + m'i)^]]. (B2) 

Note that our convention for e^^ctp^ namely, eoi23 = 1, is different from that in [lo| . 

The analytic expressions for P ^ S,A transition form factors can be obtained from that of 
P ^ P,V ones by some simple replacements. Hence, we list the explicit expressions for P ^ P 
and P ^ V transition form factors in llC 



No 
167r3 



167r3 



h' h" 

dX2Cpp'^ ^^^^ 

^X2N'iN'{ 
X2{m'i — m'l)"^ — xi{m'i 
9h' h" ( 



— ma)^ — xi{m'l — ma)^ 

xiXaM'^ — p'± — m'im2 + {m'[ — m2){x2m'i + xima) 
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+2 



q-P 



p'l + 2 



+ 2 



q \ q J q 

f'2 , o™ T\j'2 



M 



112 



X2{q^ + q-P) 



-{x2 - xi)M'^ + 2xiMo^ - 2(m'i - m2)(m'i + m'i) 



(B3) 



and 

9{q') 

f{q') 



a+iq"^ 



a-iq') 



V 



. (p± • q±f 



Nr. 



dx2d pj_ 



h' h" 

•2 / l^p'i'V 



2xx(m2 - m[){M[f + Mq^) - Axim'lM(f + 2x2^;^ • P 



167r3 J - '-X2N[N^ ^ 
+2m2q^ - 2xim2(M'2 + M"^) + 2{m[ - m2){m[ + m'l f + 8{m[ - 1712) 



+2K + + q ■ Pf-^ - 4iA±i4^ 



2xi{M'^ + Mq^) 



,2 , (p1 • 



-2(^2 + g . p) 



9^ 



2(m'^ — m'l){mi — 7112) 



dx2(fp± ^^f < (xi - X2){x2m'i + Xim2) - [2xim2 + m{ + {x2 - xi)m[] '^^ ^'^ 



167r3j - '-X2N[N'{ 

X2(^ +p'_L • g± 



X2q^Wy 



p'^ ■ A + (^i"i2 + X2m;)(xim2 - xsm'/) 



167r3 



dx2d'^p'j_ — ^T^^"! 2(2xi — 3)(x2^i^l + Xim2) — 8{m'i — 1712) 



'X2N[N^ 



P'l , ^ {PL-qA-f 

n ~r ^ A 



-[(14- 12xi)mi - 2mi 



12x,)r.2]^ 



+ J_ [ + - g2 + 2(m; - 7712 )(m;' + m2)](Af + 4'^ - 4'^) 

Wy \ 

+Z2(34'^ - 24'^ - 1) + \[xi{q^ + q-P)- 2M'2 - 2pi • 
-2m'i(m'/ + m2) - 2m2{m[ - 7712)] (vli^^ + 4^'' " 1) 



+q-P 



PX , ip'± • g^)^ 

g2 + g4 



(4^« - 3) 



(B4) 



We next give the results for p'lP'iP'iP'i and p'ip'ip'iN2- In Eq- (|3.12|) . under the typical integration 

Nr. r dX2d'^p'^ I I, 



■ h' h" 



PM 



(B5) 



167r3 J x2N'iNl 

in a P ^ M transition matrix element, p'lP'iP'iP'i in S^^'^ can be expressed in terms of three external 



momenta, P, q and a). Up to the first order in u), we have 

1=1 



3=1 k=l 



(B6) 
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where 



{gg)fj.uai3 = giiuOais + 9pLagui3 + g^ipgua, 

{gPP)fiual3 = gpivPaPlB + gfiaPuP/S + gfijsPuPa + gajsPfiPu + gufiPt^Pa + guaP^Pp, 

{gPQ)fj.ua/3 = gtMuiPaq/3 + qaPi3) + permutations, 

{gQQ)fiuai3 = gpLvQaqf} + +gfiaquqf3 + g^ipquq^ + gapq^iqu + gu/3qfj,qa + gvaqti.q(}, , 

{PPPP)fiual3 = PpiPuPaPpi 

{PPPq)^,uap = Pf^PuPaqp + Pf^PuqaPp + Pf^quPaPp + q^^PvPaPp, 
{PPqq)^uaf3 = P^iPyqaqp + permutations, 

{Pqqq)iiua(3 = P,jiquqaqp + qiiPyqaqp + qi^qvPaqp + q,iquqaPp, 
{qqq)/iuap = qf^q^qaqp, 



(B7) 



{gP^)^ival3 = 

{Pqq^)iiuap = 
{gq^)tjiuai3 = 

{qqquj)^val3 = 



Cj-P 



1 



1 



[g^_Lv{Pa^l3 + + permutations], 

(Pf^P^Padlp + Pf,PuCjaP(3 + Pf^COuPaPp + Co^PyPa^Pp), 



Cj-P 
1 



[{P^Pyqa + PfiquPa + q,iPuPa)i^p + permutations], 



tD-P 
1 



[{P/iqiyqa + q/iPuqa + q,j.qiyPa)'^p + permutations], 



Cj-P 
1 



Cj-P 



gfiuiqaojp + oJaqp) + permutations], 

{qfiquqai^p + qfiqv^aqp + qfi^uq^qp + ^^lquqaqp)- 



By contracting Pi^PipPiaP'ip with cD^, and 5"^, and comparing with the complete expressions of 
p'lfiP'ivP'ia p'l^p'ij^ shown in 0], we obtain 







4(4) 

^2 


-/l(l)/l(3) 4(4) _ 4(1) /.(3) 






— ^2 ^2 


1,(4) 


j(4) _ ,(1) ,{3) ,(4)_ 


,(1)4(3) 




— ^1 ^5 7 


^8 


_/,(l)/l{3) 4(4) _ /l(l) /l(3) 


(B8) 

q^ 




= ^S^)cf ) - 


,(4) 
^1 ) 




— ^1 -D2 ^3 ' 


^(4) 


= ^«c(^) - 


.(4) 


r-(4) _ 4(3)^(1) .Q-Pa^) 


r^(4) _ 4(3)^(1) ,o5--P 4(4) 
"-^2 ~ ^6 ^1 ~'~ ^2 4 ' 


where (lo| 











^ 2 ' 



A. 



(1) 



^(1) _ Pi_ -qi. ^ ^(1) 



2 ^2 



-iV2 + ^2, 



p± • qi. 



A 



Z2 = N[ + m'{ -ml + {l-2xi)M'^ + {q^ + q-P) 
iP± ■ ^(2) 



(2) 



/2 
-P± 



iA\ 



(1)a2 



A 



(2) 



(2) 



(A 



1,(2) 



,(3) _ 4(1) ,(2) 



(3) 



,(!),(!) 
^1 ^2 ' 



,(1) ,(2) 
^2 ^1 ' 



(B9) 
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zl (3) - 4(1)^(2) 4(3) _ /.(I) 4(2) 4(3) _ 4(1)4(2) 
^3 — ^1 ^2 ' ^4 ~ ^2 ^2 ' ^5 ~ ^1 ^4 ' 

4(3) _ 4(1)4(2) 1^(1)4(2) 

yig — 74.2 g2 2 ^1 ' 

Following the prescription in ^Ifl], the spurious contributions cj^2 should be vanished by including 
the zero mode contribution and we have 

Af'>N,^A^?Z,+ ^-^A^^\ AfN,^Afz, + ?,'l^Af, (BIO) 

which lead to the p'^p'i,p'aN2 formula shown in Eq. ()3.14p . Note that in general are non- vanishing 
by themselves, but they do vanish under integration in some choice of vertex function There 
are some attempts to include these effects for generic vertex functions 14]. The important of these 
effects can be checked numerically. For example, we have checked numerically that the integral of 
Eq. (|B5|) with S^^'^ replaced by b\^'^ are vanishingly small. In practice, one only needs terms 
for p'l - ■ ■ p'l formulas. 
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